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short  time  tracking  of  planets,  if  so  desired.  The  details  of  this  program 
were  presented  in  the  interim  report,  (PAR  Report  #75—32)  submitted  in 
December  1975. 

^There  are  two  separate  systems  for  tracking  satellites.  In  the  absence 
of  any  satellite  differential  position  measuring  equipment,  Simplified 
Generalized  Perturbation  theory  is  used  to  predict  satellite  position. 

Using  NORAD  satellite  card  parameters  and  a timing  input,  the  program  is 
sequenced  to  compute  the  satellite's  direction  and  outputs  these  angles  to 
the  mount  servo  systems.  In  the  second  system  a computer-TV  interface  was 
fabricated  to  measure  the  satellite  differential  position  after  obtaining  it 
in  the  field  of  view.  This  measured  position  is  fed  back  to  the  computer, 
which,  in  turn,  uses  that  as  observed  error  data  for  a Kalman  filter.  The 
filter  calculates  updated  parameters  while  tracking  the  satellite  position 
as  a function  of  time. 

The  report  includes  a detailed  technical  description,  program  documentation, 
and  a system  user's  manual.  __ 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


1.1.  INTRODUCTION 

This  report  describes  the  technical  aspects  of  software  developed  by 
Pattern  Analysis  and  Recognition  Corporation  (PAR)  for  driving  the  optical 
tracking  mount,  at  RADC's  AOTF  (Advanced  Optical  Test  Facility),  Verona  test 
annex  with  a NOVA  800  computer.  The  software  developed  to  drive  the  mount 
for  tracking  stars  was  described  in  PAR  Report  No.  75-32  (Contract  */F30602- 
75-C-0144);  the  hardware  designed  to  interface  the  computer  with  the  mount 
was  described  in  another  report  titled  "Coelostat  Tracking  Mount/NOVA  800 
Hardware  Digital  Interface"  (PAR  Report  No.  75-29).  This  report  describes 
the  technical  aspects  of  the  software  developed  to  track  satellites  of  known 
orbital  parameters. 

Passive  satellite  tracking  routine  makes  use  of  Simplified  General 
Perturbation  theory  in  conjunction  with  recently  computed  orbital  parameters 
to  predict  the  position  of  the  satellite  at  any  given  instant  of  time.  The 
section  titled  "ACQUISITION"  describes  a spiral  search  routine  designed  to 
perform  a systematic  search  around  any  calculated  satellite  position;  this 
will  be  of  great  use  when  inexactness  of  satellite  parameters  causes  the 
mount  to  point  away  from  the  true  direction  and  it  becomes  necessary  to 
acquire  the  satellite  in  the  field  of  view.  Active  satellite  tracking 
describes  the  software  designed  not  only  to  predict  satellite  position,  but 
also  to  apply  suitable  corrections  to  its  orbital  parameters  according  to 
the  observed  azimuth-elevation  angle  information  of  the  satellite.  The 
routine  uses  estimation  techniques  in  the  form  of  a linearized  Kalman  Filter. 
The  hardware  designed  to' measure  the  differential  azimuth  - elevation  angles 
is  described  in  a later  section  of  this  report.  Program  descriptions  and  a 
User's  Manual  are  also  included. 

1.2.  THE  TRACKING  MOUNT 

The  Coelostat  tracking  mount  permits  an  observing  instrumentation  system 
to  remain  stationary  while  the  tracker  points  anywhere  within  the  hemisphere. 
This  particular  system  allows  excursions  of  approximately  ten  degrees  below 
the  horizon  with  one  meter  clear  aperture  (Figure  1-1). 


The  mount  is 
that  is  rotatable 
azimuth  rotationa 
the  line  of  sight 
tion  of  the  mount 
the  azimuth  axis 
mately  200°.  The 
azimuth  axis  and 


a two-mirror  system.  Both  mirrors  are  mounted  on  a platform 
to  any  azimuth  position.  One  mirror  is  centered  on  the 
l axis  and  stationary  with  respect  to  the  platform;  it  folds 
ninety  degrees  to  become  perpendicular  to  the  axis  of  rota- 
. The  second  mirror,  the  elevation  mirror,  is  displaced  from 
and  cradled  to  rotate  through  the  elevation  range  of  approxi- 
Coelostat  is  driven  by  two  independent  servos , one  for  the 
the  other  for  elevation  (Figure  1-2).  By  selecting  "external 
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position  designate"  mode,  the  mount  can  be  driven  by  the  digital  computer. 

To  orient  the  line  of  sight  in  a desired  direction,  the  software  supplies 
the  corresponding  azimuth  and  elevation  angles  to  the  mount  through  suitable 
interface  hardware.  For  detailed  description  of  the  mount  refer  to  Refer- 
ence [1]. 


1.3.  CAPABILITIES  OF  THE  SOFTWARE 

The  STARTRACK  program  requires  the  current  right  ascension  and  declina- 
tion of  the  star  to  be  tracked,  year  and  day  number  as  operator  inputs.  The 
clock  time  input  has  been  hardwired.  Stars  above  an  elevation  of  approxi- 
mately 15°  will  almost  always  be  visible  in  the  center  of  the  field  of  view 
of  the  main  telescope,  weather  conditions  permitting.  Up  to  seven-magnitude 
stars  have  been  photographed.  Visibility  of  stars  at  low  elevations  cannot 
be  guaranteed  because  of  haze  and  atmospheric  refraction.  The  star  catalog 
[2]  lists  the  RA  and  DEC  values  of  stars  brighter  than  seven-magnitude. 

Since  the  catalog  tabulates  the  apparent  traces  of  the  stars  at  each  ten-day 
interval,  it  is  unnecessary  to  account  for  precession  and  nutation  of  the 
ecliptic  in  the  program.  These  phenomena  have  been  considered  while  tabu- 
lating the  catalog. 

The  SATTRACK  program  requires  the  NORAD  two  card  satellite  parameter 
set  [3]  as  input  along  with  the  day  and  the  year.  If  the  parameter  set  has 
been  recalculated  recently  (epoch  time  within  the  past  week  approximately), 
the  satellite  will  often  be  visible  when  it  is  out  of  earth's  shadow,  at 
least  in  the  viewing  telescope.  The  viewing  telescope  has  its  optical  axis 
parallel  to  that  of  the  main  telescope  and  looks  at  the  sky  through  the 
coelostat.  Because  of  its  larger  field  of  view  (one  degree  as  opposed  to 
five  arc-minutes  of  the  main  telescope),  satellites  not  in  the  main  scope 
may  be  visible  in  the  viewing  scope.  The  satellite  can  then  be  brought  to 
the  center  of  the  field  of  view  using  the  Mount  Control  routines  described 
in  the  User's  Manual.  The  program  is  a simplified  version  [3]  of  the  perturba- 
tion theory  of  Kozai  [4]  and  it  does  not  compute  any  corrections  attributed 
to  any  of  the  following  causes. 

1.  Atmospheric  drag 

2.  Solar  pressure 

3.  Lunar  and  Solar  gravitation  effects 

4.  Ray  bending  in  the  atmosphere 

5.  Structural  deformations  of  the  telescope  and  mount 

6.  Alignment  errors 

A number  of  satellites  were  tracked  successfully  with  the  passive 
tracking  program.  While  tracking  a satellite,  the  tracking  rate  never 
exceeded  the  slew  rate  of  the  mount. 
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SECTION  2 


USER'S  MANUAL 


2.1.  INTRODUCTION 

This  user's  manual  is  written  with  the  inexperienced  user  in  mind.  All 
commands  have  been  grouped  according  to  their  functions  and  their  explanations 
have  been  augmented  by  examples.  Unless  otherwise  specified,  all  the  commands 
described  in  this  manual  perform  identical  operations  in  SATTRACK  and  STAR- 
TRACK  programs.  It  is  assumed  that  an  experienced  mount-operator  will  be 
responsible  for  manning  the  mount-console. 

All  commands,  except  those  which  are  at  the  RDOS  system  level  (RDOS, 
SATTRACK,  STARTRACK  and  CARDINPUT),  are  comprised  of  four  letters.  The 
system  will  accept  just  four  letters,  but  the  description  of  commands  includes 
more  than  just  four  letters  so  as  to  enable  the  reader  to  associate  the 
command  with  the  function  it  performs. 

An  alphabetical  index  of  commands  is  included  in  Section  2.6  for  cross- 
reference  . 


2.2.  PRELIMINARY  OPERATIONS 

2.2.1.  System  Initialization 

The  program  is  controlled  by  the  Real-time  Disk  Operating  System  (RDOS) 
on  the  Nova  800  computer  system  and  the  following  initialization  procedure 
must  be  adopted. 


2. 2. 1.1. 


Insert  the  disk  with  the  required  routine  (STARTRACK  and  SATTRACK) 
into  the  disk-drive  and  turn  "on"  the  power  switch  key  located  above  the 
disk-drive. 


2. 2. 1.2. 


Press  the  RUN-LOAD  switch  to  "RUN"  and  wait  till  the  READY  light  comes 
on. 


2. 2. 1.3. 


Turn  on  the  power  to  the  CPU, TTY  and  the  line  printer,  and  also  the 
'SELECT'  button  on  the  line  printer.  The  power  for  the  card  reader  is 
located  on  the  back  of  the  card  reader. 


2. 2. 1.4. 


Load  the  octal  number  100033  (1000000000011011  in  Binary)  in  the  toggle 
switch  register  on  the  front  of  the  CPU  and  press  program  'LOAD'  up. 
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2. 2. 1.5. 


The  TTY  should  print  'FILENAME?1  the  proper  response  to  which  is  RDOS 
followed  by  a carriage  return  ( J ) . 

Operator  commands  or  responses  are  generally  followed  by  a CR  ( J ) 
and  system  responses  have  been  underlined  to  avoid  confusion. 


2. 2. 1.6. 


If  the  system  responds  with  ’PARTITION  OCCUPIED.  TYPE  C TO  CONTINUE' , 
then  type  only  'C'  not  followed  by  carriage  return. 


2. 2. 1.7. 

Answer  the  questions  regarding  date  and  time  by  RDOS,  which  will  then 
respond  with  'R_'  on  the  TTY.  The  'R_'  indicates  that  RDOS  is  in  control  and 
R_EADY  to  accept  commands.  This  completes  the  sequence  of  operations  for 
system  initialization. 

2. 2.1.8. 

To  leave  RDOS  and  shut  off  the  system,  type  l?0UT@  & . The  system 
responds  with  'MASTER  DEVICE  RELEASED.'  Then  turn  off  the  power  on  all  the 
devices. 

2.2.2.  Track  Initialization 

2. 2. 2.1.  PRINTMANUAL 

In  case  a printout  of  the  tracking  manual  describing  all  the  commands 
is  required,  type  PRINTMANUAL  ii . This  command  prints  the  complete  manual 
on  the  line  printer. 

2. 2. 2. 2.  CARDINPUT 

To  update  the  satellite  list,  type  CARDINPUT  i and  answer  the  questions 
and  read  iu  the  cards. 

2 . 2 . 2 . 3 . STARTRACK 

To  track  stars  or  any  other  heavenly  body  with  a slowly  varying  right 
ascension  and  declination,  type  STARTRACK • 

2. 2. 2.4.  SATTRACK 

To  initiate  satellite  tracking,  type  SATTRACK  1 1 and  then  GETP J to  get 
the  parameter  of  the  satellite  to  be  tracked  (Section  2. 3. 2.1.).  Example:# 

I When  the  system  prints  an  asterisk  it  means  that  the  SATTRACK  main 
command  line  interpreter  (CLI)  is  in  control,  as  opposed  to  R^ which  indicates 
control  by  RDOS  CLI. 
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. . X 


R 

SATTRACKi 

INPUT  THE  DAY  AND  YEAR  (13,  IX,  12) 

325,76 

A, 

GETPjA 

LIST  NUMBER 

3* 

cont^ 

2. 2. 2. 4.1.  KALCARD 

R 

KALCARD^.  This  is  a command  executable  at  the  RDOS  level  only. 

This  routine  allows  you  to  store  new  cards  to  the  disk  in  a 
file.  Example: 

R 

KALCARD  * 

SATELLITES  MUST  BE  DELETED  IN  THE 

ORDER  IN  WHICH  THEY  APPEAR  IN  THE  LIST. 

FIRST  SATELLITE  TO  BE  DELETED.  The  system  asks  the  operator  to  delete  a 
satellite  if  he  wishes.  If  the  operator  answers  with  a $ , no 
satellites  will  be  deleted. 

4^  The  operator  requests  satellite  #4  to  be  deleted. 

NEXT  SATELLITE?  The  system  asks  for  the  next  satellite  to  be  deleted. 

Oil  The  operator  specifies  none 

INPUT  THE  NUMBER  OF  CARDS  TO  BE  READ.  The  system  asks  for  a maximum  number 
of  cards  to  be  read. 

lOOOOjt  The  operator  specifies  all  of  them. 

The  routine  will  keep  reading  cards  until  10000  cards  have 
been  read  or  until  an  end  of  file  card  has  been  reached  (not 
supplied  by  NORAD).  If  no  end  of  file  card  is  available,  turn 
off  the  card  reader. 

DO  YOU  WISH  THE  WORST  CASE 

PARAMETERS  TO  BE  KEPT?  The  system  asks  the  operator  if  he  wants  the  Kalman 
filter  parameters  to  be  the  worst  case  of  the  past  history  of 
the  satellite  to  be  the  ones  kept  for  tracking.  The  usual 
response  is  'YES’.  After  many  updates,  the  answer  should  be 
no  for  one  update  then  yes  again. 

YES*1  The  operator  specifies  a yes. 

the  question 

"INPUT  THE  DIAGONAL  TERMS  FOR  THE  PK  AND  QX  MATRIX." 

The  operator  should  then  give  all  (10  5)'s  if  he  knows  nothing 
about  the  filter.  A more  intelligent  guess  of  the  parameters 
would  be  a satellite  of  nearly  the  same  card  set  and  use  the 
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diagonal  terms  for  that  satellite.  These  terms  are  the  second 
set  of  parameters  printed  on  the  line  printer  for  a satellite. 
The  twelve  elements  could  be  entered  all  on  one  line  or  on  any 
number  of  lines.  The  reason  for  this  question  is  that  this  is 
the  first  known  card  set  for  this  satellite,  since  no  past 
history  is  known  for  it.  The  operator  is  then  given  a chance 
to  guess  at  it. 

LOAD  $CDR?  TYPE  ANY  KEY  WHEN  READY  The  operator  is  to  load  the  card 
reader  with  the  NORAD  cards. 

^ The  operator  is  ready  to  read  in  the  cards. 

MORE  CARDS?  The  card  reader  has  reached  an  end  of  file  or  has  been  turned 

off.  The  operator  can  read  in  more  cards  by  typing  "YES". 

NOv.  He  has  none. 


INPUT  AN  8 CHARACTER  FILE  NAME  FOR  AN 

INPUT  FILE.  The  system  wants  to  know  where  to  store  the  newly  created  list 
of  parameters.  This  file  must  already  be  present  on  the  disk 
and  must  be  writable. 


KALLIST ^ The  operator  specifies  the  file  "KALLIST"  as  the  output  file. 

The  trailing  spaces  are  to  pad  the  word.  "KALLIST"  is  the 
list  of  parameters  used  for  the  Kalman  filter. 


STOP 

R_  The  operator  is  back  at  RDOS. 

2. 2. 2. 4. 2.  KALTRACK 


Main  Kalman  filter  tracking  routines. 

2. 2. 2. 4. 3.  FILTER 

FILTER^  This  command  implements  a Kalman  filter  on  the  present  set  of 

card  parameters  being  used  for  tracking.  The  output  of  this 
command  is  a new  set  of  card  parameters.  The  data  resulting 
from  the  filter  is  printed  on  the  line  printer.  The  operator 
is  asked  the  question,  "AZ,EL  OBSERVED  ERRORS?"  The  operator 
is  being  asked  for  any  offsets  he  wants  added  on  top  of  the 
offsets  already  present  in  the  system.  The  system  then 
responds  with  the  question:  "DO  YOU  WISH  TO  UPDATE  THE  PK 

MATRIX?"  The  PK  matrix  is  the  system  confidence  matrix!  The 
normal  answer  to  this  question  is  "NO".  An  understanding  of 
the  filter  would  be  needed  to  answer  yes.  Example: 

A 

FILTER  d 

AZ , EL  OBSERVED  ERRORS ■ (D ,M ,S ) 

DO  YOU  WISH  TO  UPDATE  THE  PK  MATRIX? 


NOii 


* The  Kalman  filter  has  been  executed,  the  system  t^me  offsets 

and  mount  offsets  have  been  used  as  the  error  inputs  to  the 
filter,  and  the  offsets  have  been  cleared  from  the  system. 

2.2.2. 5.  Q TRACK 

QTRACKjZ  Queue  the  system  to  start  and  stop  tracking  at  operator- 

specified  time  instants  in  universal  time  (U.T.).  Example: 

*_  Command  Line  Interpreter  (CLI)  is  in  control 

QTRACKi  Operator  wishes  to  queue  the  system 

START  TIME,  FINISH  TIME  (UT,  H,  M,  S)  The  system  asks  for  start  and  end 
times. 

15,5,8 J Operator  wants  the  track  to  start  at  15h.  5m.  8sec.  U.T. 

15,11,8  iZ  The  track  is  to  stop  at  15h,  11m,  8sec,  U.T. 

AZ,  EL  START  POSITIONS  = XXX.XXXXX  XX . XXXXXX 


AZ,  EL  FINISH 


specified  start  time 
POSITIONS  = XXX.XXXXX 


XX. XXXXXX 


specified  end  time 

Control  has  been  relinquished  to  CLI 


Object  position  at 
Object  position  at 


the 


2. 2. 2. 6. 


QUIT 


QUITjZ 


This  command  stops  the  track  and  returns  control  to  RBOS. 
If  newly  added  parameters  to  the  list  are  to  be  saved,  the 
list  must  be  written  on  the  disk  before  issuing  this  command 
in  STARTRACK  or  SATTRACK.  KALTRACK  is  always  operating  from 
the  disk.  Example: 


QUITJ 


STOP 

R 


Operator  wishes  to  stop  the  track 
The  command  is  executed 

RDOS  indicates  that  the  tracking  routine  is  closed 
RDOS  is  in  control 


2. 2. 2. 7.  START 


STARTS 


Start  the  tracking  routine  again.  This  has  the  same 
effect  as  coming  in  from  RDOS  except  that  the  parameter  list 
contains  any  additions  made  since  the  last  disk  WRITE  . All 
flags  have  to  be  reset  by  the  operator.  Example: 

CLI 


STARTS 

INPUT  THE  DAY  AND  YEAR  (13,  IX,  12) 
166,761  Day  #166  in  the  year  1976 

* CLI 
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2. 2. 2. 8.  STOP 


STOP  > 


Track  is  stopped  but  control  is  retained  by  the  tracking 
program.  This  allows  the  operator  to  start  a new  track  or 
have  the  system  waiting  for  him. 


2. 2. 2. 9.  CONT 

CONT  This  command  causes  the  system  to  continue  tracking  by 

computing  new  AZ,EL  angles  and  transferring  them  to  the 
mount.  This  is  an  active  tracking  mode  and  any  searches  can 
now  be  performed. 

2.3.  DATA  STORAGE  AND  RETRIEVAL 

Data  required  for  tracking  is  stored  in  a star  or  satellite  list  which 
can  be  updated  as  necessary. 


2.3.1. 


Storage 


2. 3. 1.1.  XFER 


XFER  1 


line  1: 


line  2: 


Transfer  the  given  set  of  parameters  to  the  appropriate  list. 

The  data  format  is  different  for  a star  list  and  a satellite 
list.  In  a star  list  the  parameters  are  star  #,  brightness, 
right  ascension  (hours,  minutes,  seconds),  declination  (degrees, 
minutes,  seconds),  so  that  the  units  of  RA  and  DEC  are  the 
same  as  those  in  the  star  catalog. 

The  parameters  for  satellite  list  are  in  the  same  order  and 
units  as  on  the  cards  ('two  card  format')  and  must  be  entered 
on  two  lines.  Order  of  parameters: 

Epoch,  year,  revolution  #,  mean  motion  (n),  inclination  (i), 
eccentricity  (e) 

Mean  anomaly  (M),  RA  of  the  ascending  mode  ( 0 },  argument  of 
perigee  (w),  first  derivative  of  mean  motion  (n),  second 
derivative  of  mean  motion  (n),  satellite  #.  This  command  is 
not  valid  in  KALTRACK. 


2. 3. 1.2.  READ 


Delete  the  list  in  core  and  replace  it  with  the  list  on  the 
disk.  (Not  valid  in  KALTRACK.) 


2. 3.1. 3.  WRITE 


WRITE  )■  Write  the  list  in  core  over  the  list  on  the  disk.  The  list 

in  core  is  left  unmodified.  No  copy  of  the  overwritten  list 
on  the  disk  is  retained.  (Not  valid  in  KALTRACK.) 


2. 3.1.4.  STORE 


STORED  This  command  stores  the  filtered  card  parameters  out  to  the 

disk  list  called  "KALLIST".  The  system  asks  for  "LIST#".  The 
system  stores  the  card  parameters  in  the  location  specified 
unless  the  specified  number  is  greater  than  the  final  number 
in  the  list.  In  this  case,  the  parameters  are  stored  in  the 
slot  after  the  last  card  satellite  in  the  list.  Example: 


The  newly  created  card  parameters  are  stored  in  the  last 
vacant  position  in  the  list  because  the  list  does  not  contain 
1000  satellites. 

2.3.2.  Retrieval 

2. 3. 2.1.  GETP 

GETPii  Get  parameters  from  the  parameter  list  and  write  them  on  the 

line  printer  and  the  system  goes  in  a halt  mode.  Tracking  can 
now  be  started  by  the  CONT  ^ command.  Example: 

GETP4 

LIST  NUMBER?  (12)  The  system  wants  to  know  the  place  in  the  list  of 
parameters. 

3^  The  operator  wishes  to  track  the  object, 

(star  or  satellite  depending  on  the  program  in  control)  #3  in 
the  list. 

Note:  This  number  is  not  the  same  as  the  star  catalog  number 

or  the  satellite  number  but  only  the  number  in  order,  in  the 
list. 

2. 3. 2. 2.  LIST 

LIST4  This  command  prints  out  the  list  resident  in  the  core  to  the 

!*  terminal  without  any  more  responses  from  the  operator.  (Not 

S'  • valid  in  KALTRACK. ) 

I ' 

2. 3. 2. 3.  PRINT 
PRINT 4 


V 

5 


Print  a listing  of  the  mission  to  the  line  printer.  For  star 
track  the  print  reads  as  follows:  year,  day,  hour,  minute, 

seconds,  output  azimuth,  output  elevation,  azimuth  read  from 
the  encoders,  elevation  read  from  the  encoders,  and  number  o£ 
outputs  since  the  last  print  to  the  line  printer. 


* 

ST0RE4 

LIST#. 

10004 


For  satellite  track,  the  printout  is:  year,  day,  hour, 

minute,  seconds,  true  anomaly,  output  azimuth,  output  eleva- 
tion, azimuth  read  from  encoders,  elevation  read  from  encoders, 
number  of  outputs  since  the  last  print  to  the  LP,  and  number 
of  clock  errors  since  the  last  format  to  the  LP. 

The  system  also  allows  the  operator  to  choose  the  print  mode 
he  wants,  the  print  interval  in  seconds,  and  the  number  of 
prints  per  camera  exposure.  Example: 


PRINT 

INPUT  PRINT  MODE  System  wants  the  print  mode. 

1 = NO  PRINT  If  any  printing  is  being  done,  stop  it. 

2 - DISK  PRINT  Write  the  data  on  the  disk. 

3 = $LPT  PRINT  Print  the  data  on  the  line  printer. 

4 = PRINT  BOTH  Write  on  the  disk  and  the  LPT. 

3 ii  Operator  chooses  to  print  on  the  LPT. 

PRINT  INTERVAL 

10*.  Data  is  to  be  printed  every  ten  seconds. 

PRINT  INTERVAL  PER  CAMERA  QUEUE 

10000*!  This  is  a large  number  which  indicates  that  no  pictures  are  to 

be  taken. 

* CLI 

2. 3. 2.4.  HEADER 


HEADER i Output  a header  to  the  line  printer  just  before  the  next  print 

which  must  be  less  than  or  equal  to  eighty  characters. 

Example : 


HEADER*: 

HEADER'. 
THIS  IS  A 
ft 


Operator  wishes  to  insert  a header  in  the  printout  on  the  line 
printer . 

System  is  ready  to  accept  the  message. 

TEST.l  Message  to  be  inserted  in  the  printout. 

Message  is  inserted  and  control  returned  to  CLI.  System 
offsets  of  time  and  mount  position  are  printed  also. 


2.4. 


MOUNT  CONTROL 


This  section  describes  the  commands  used  to  direct  the  mount  to  desired 
orientations,  either  under  operator  control  or  under  a track-routine  control. 

2.4.1.  Operator-Specified  Mount  Control 

2. 4. 1.1.  MMOUNT 


MMOUNTj) 


Manually  position  the  mount  by  using  the  TTY  to  add  or  sub- 
tract a preset  set  of  increments  to  the  mount  offsets  already 
in  the  system.  If  the  mount  is  still  under  control  of  the 
command  CMOUNT,  it  is  overridden.  There  are  subcommands 
available  for  the  manual  control.  (Not  valid  in  KALTRACK.) 
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Av 

E». 

>2  or  . jt 
<2  or 
?v 

i 


Control  the  azimuth  axis 
Control  the  elevation  axis 

Add  a positive  offset  to  the  axis  under  control. 

Add  a negative  offset  to  the  axis  under  control. 

Leave  this  mode  and  give  the  total  offsets  in  degrees. 
Leave  this  mode  without  giving  the  offsets. 


Any  other  character  that  is  typed,  is  ignored.  Example: 


MMOUNT 

AZ,  EL  SLIP  INC.  (D,M,S)  PER  < OR  > System  axis  for  increment  (decrement) 


value  per  increment  (decrement)  command 

3,30,0^ 

Azimuth  increment  is  3.5  degrees. 

5,15,0(2 

Elevation  increment  is  5.25  degrees 

E J 

Decrement/increment  the  elevation  axis. 

>d 

Add  one  positive  increment. 

. k 

Add  one  more  increment. 

• k i 

Add  one  more  increment . 

Subtract  one  increment. 

<i 

Subtract  one  more  increment. 

Ad 

Decrement/increment  the  azimuth  axis. 

>« 

Add  one  offset 

< i i 

Subtract  one  offset. 

, d 

Subtract  one  more. 

< d 

Subtract  one  more. 

< ^ 

Subtract  one  more. 

. k. 

Add  one  offset. 

AZ,  EL  OFFSETS 

IN  DEGREES  ARE  - XX.XXXX  XX.XXXX« 

» cli 

2.4.1. 2.  C MOUNT 


C MOUNT 


• i or  0 ) 


>n‘ 

. u 
U 
d 

id 


This  command  achieves  continuous  mount  positioning  according 
to  the  increment  rate  per  second  supplied  by  the  operator. 

Same  subcommands  as  those  under  MMOUNT  are  available  here  also 
to  facilitate  tracking,  but  the  increment  (decrement)  subcom- 
mands achieve  slightly  different  results  under  this  continuous 
mount  positioning  mode. 

Reverse  the  direction  of  increment  only  for  the  axis  under 
control. 

Slow  down  the 
Slow  down  the 
rate. 

Speed  up  the  axis  under  control  to  twice  the  current  rate. 
Speed  up  the  axis  under  control  to  thrice  the  current  rate. 
Stop  the  axis  under  control. 

System  continues  the  operation  but  returns  to  CLI. 

Leave  this  mode,  print  to  the  operator  the  total  system 
_2ff.sets  and  stop  any  CMOUNT  positioning  activity. 


axis 

axis 


under  control 
under  control 


to  half  the  current  rate, 
to  a third  of  the  current 


ft  These  reported  values  are  of  total  offsets  and  not  only  those  due  to 
the  current  MMOUNT  command,  if  some  offsets  already  existed  before  ini- 
tiating MMOUNT. 


2-9 


»»  ■**•  •**»  • » * \ 


v • K7- 


If  CLI  is  under  control,  and  current  mount  positioning  is  to  be  altered, 
this  mode  is  re-entered  by  typing  CMOUNT  d which  causes  the  axes  movements  to 
stop.  If  any  character  other  than  the  above  is  typed,#  then  the  system 
continues  to  move  the  mount  at  the  latest  chosen  rate.  Example: 

* 

CMOUNT^ 

AZ,  EL  SLIP  RATE  (D,M, S/SEC)?  The  system  asks  for  the  rate  in  (degrees, 

minutes,  seconds) /sec.  of  time. 

Azimuth  rate  would  be  2°/sec . 

Elevation  rate  would  be  l°/sec. 

Position  the  elevation. 

Start  the  positioning  at  the  given  rate. 

Double  the  elevation  changing  speed. 

Reduce  the  elevation  speed  to  a third. 

Reverse  direction  of  motion. 

Stop  elevation  positioning. 

Start  azimuth  positioning  at  the  given  rate. 

Slow  down  to  a third  of  the  speed. 

Stop  azimuth  change,  give  offsets,  and  return  to  CLI. 

AZ,  EL  OFFSETS  IN  DEGREES  ARE  - XX.XXXXXXX  XX.XXXXX 

...  ~ - - — — - — 1 ------  — — 

2. 4. 1.3.  INCM 

INCM  ' Increment  the  previous  mount  offsets  by  the  given  amount. 

(Not  valid  in  KALTRACK.)  Example: 

* 

Tncm^ 

AZ,  EL  - INC  (D,M,S ) 

2.3.4,  Azimuth  increment  is  2°,  3 minutes,  and  4 seconds. 

4,5,6/  Elevation  increment  is  4°,  5 minutes,  and  6 seconds. 

AZ  OFFSET  = X.XXXX  EL  OFFSET  = X.XXXX 

7: 

2. 4. 1.4.  SETM 

SETM  Set  the  mount  offsets  to  operator-specified  new  offsets  by 

deleting  previous  offsets.  Example: 

ft 

SETM  *■ 

AZ,  EL  (D,M,S) 

3,2,4; 

5,6,7. 

ft 

# In  the  example,  since  L is  "any  character  other  than  above,"  it  has  the 
effect  of  starting  the  positioning  at  the  given  rate. 


2 .0 . 01 

1.0.  01 
Ev 

l a 


> y. 


0* 
/* 
A * 
L J- 
A 

» r 

•5  > 
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2.4.2. 


Track  Routine  Mount  Control 


In  the  previous  subsection,  mount  control  was  achieved  by  mount  offsets 
directly  specified  by  the  operator.  The  commands  did  not  go  through  any 
tracking  computations.  Both  commands  in  this  subsection  are  valid  for  star 
track  routine  only.  They  may  be  used  to  center  a star  in  the  field  of  view 
by  altering  the  RA  and  DEC  values  instead  of  AZ  and  EL  angles. 

2. 4. 2.1.  INCS 

INCS  This  is  a valid  command  in  the  star  track  program  alone  and 

it  permits  the  operator  to  add  or  subtract  an  increment  to  the 
RA  and  DEC  of  the  star  already  being  tracked.  The  increment 
supplied  does  not  alter  the  parameters  in  the  list  which  can 
be  changed  only  by  the  'XFER'  command.  Example: 

* 

Incs  > 

RA,  DEC  - INC  ( H , M , S D,M,S) 

4,6,8* 

6,8,9^ 

RA  OFFSET  = XXX.XXXX  DEC  OFFSET  X.XXXX 

x — ' ■ ' ' — • — ■ ■ 

2.4. 2. 2. 

SETS  J 


SETS  ^ 

RA,  DEC  (H,M,S  D,M,S) 

1,2,3* 

4,5,6* 


SETS 

This  command  also  works  for  STARTRACK  only,  and  it  allows  the 
operator  to  specify  completely  new  RA  and  DEC  values  for  the 
star  being  tracked  without  changing  them  in  the  parameter 
list.  The  mount  starts  tracking  according  to  the  new  para- 
meters. Example: 


2.5.  SEARCH  COMMANDS 

Strictly  speaking,  a search  may  be  conducted  by  using  any  of  the  appro- 
priate mount  control  commands  described  in  the  preceding  subsection.  However, 
the  search  process  will  then  involve  a tedious  sequence  of  increments  speci- 
fied by  the  operator,  and  it  is  not  always  possible  for' the  operator  to 
select  a correct  set  of  increments.  Especially  for  the  case  when  the  object 
is  not  visible  in  the  field  of  view,  a spiral  search  can  be  conducted  by  the 
commands  described  in  2.5.1.  By  experience,  it  has  been  determined  that 
satellite  tracking  errors  generally  involve  time  offsets.  That  is,  the 
calculated  orientation  of  the  mount  may  point  slightly  ahead  or  behind  the 
actual  satellite  position  in  its  orbit.  If  such  is  the  case,  the  time  slip 
commands  described  in  subsection  2.5.2.  will  be  invaluable. 
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2.5.1. 


Spiral  Search  (Not  Valid  in  KALTRACK) 


To  use  this  option  efficiently,  the  operator  must  have  certain  familiar- 
ity with  the  mechanism  of  search.  A detailed  description  of  the  options  is 
available  in  the  section  on  acquisition  and  will  not  be  repeated  here. 


2. 5. 1.1.  SEARCH 


SEARCH  i. 


This  command  initiates  the  procedure  to  implement  a spiral 
search  in  the  line  of  sight  coordinates  with  approximately 
constant  linear  speed.  The  system  asks  a number  of  questions 
to  get  a specification  of  the  spiral  to  be  drawn  out  about  the 
calculated  azimuth-elevation  direction.  In  response  to  the 
questions,  the  operator  must  specify  the  view  time,  percentage 
of  field  of  view  overlap,  and  maximum  angle  of  search. 

Example : 


SEARCH  / 

HELP 
0/ 

INPUT  VIEW  TIME 


The  system  asks  if  the  operator  needs  help. 

The  operator  declines.  (No  other  option  is  yet  operational.) 
IN  SECONDS  This  is  the  time  for  which  a point  on  the 


spiral  is  required  to  stay  in  the  field  of  view. 

10/  View  time  is  10  seconds. 

PERCENTAGE  OF  EIELD  OF  VIEW  TO  BE  USED  This  is  the  percentage  of  the  aperture 
by  which  the  spiral  expands  after  going  through  360°. 

20/  Operator  specifies  20%  or  1 minute  since  the  angular  field  of 

view  is  5 minutes. 

MAXIMUM  RADIUS  TO  BE  SEARCHED 

5>i  Radius  is  specified  in  terms  of  angle  in  5 arc  minutes. 

TIME  TO  COMPLETE  THE  SPIRAL  IS  (IN  SECONDS)  157.5762 

? The  $ sign  indicates  that  the  search  CLI  is  in  control  instead 

of  the  main  CLI,  and  the  system  will  now  be  able  to  accept 
commands  valid  in  the  search  routine.  The  spiral  has  now  been 
initiated  and  the  mount  will  continue  to  rotate  about  the 
continuously  varying  calculated  position  till  it  reaches  the 
maximum  specified  limit  or  is  interrupted  by  the  commands  to 
be  described  presently.  The  system  does  not  respond  with  ^ 
since  the  main  CLI  is  no  longer  in  control  and  the  search  CLI 
is. 


2. 5. 1.2. 
PAUSE  J. 


PAUSE 


This  command  halts  the  radius  and  angle  changing  at  the 
instant  the  command  was  given  and  the  mount  continues  tracking 
with  the  current  offsets. 


2. 5.1. 3. 
BACK  A 


BACKTRACK 


This  command  allows  the  exact  same  spiral  to  be  traced  backward 
at  an  operator-specified  rate,  after  the  "PAUSE"  was  made  in 
the  search.  The  operator  is  asked  to  specify  the  backtrack 
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rate  as  a percentage  of  the  rate,  before  PAUSE  was  issued.  If 
no  pause  was  issued,  one  is  automatically  inserted  in  the 
command  stream.  Example: 


The  system  requires  the  backtrack  rate. 

The  operator  wants  the  system  to  describe  the  spiral  backwards 
at  half  the  previous  rate. 

If  the  response  was  -100  A instead  of  50^-  , the  spiral  would 
have  continued  to  be  described  in  the  forward  direction  at 
the  same  rate  as  was  before  the  PAUSE  command.  If  a backtrack 
is  allowed  to  continue  to  the  center  of  the  spiral,  the  search 
is  automatically  paused. 

2. 5. 1.4.  EXIT 

EXIT.A  This  command  causes  the  system  to  exit  the  search,  leaving  the 

current  offsets  in  the  system  and  printing  them  out  to  the 
terminal  in  degrees. 

2 . 5 . 1 . 5 . ABORT 

ABORT )•  The  search  is  aborted  and  any  offsets  that  were  introduced 

during  the  search  are  nullified.  Any  offsets  that  existed 
before  the  search  are  left  unchanged. 

2. 5.1.6.  CSEARCH 


$ 

BACK  j) 
RATE 
50 
$ 


I 


* 


* 

%» 

4 

\ 

t 

J 

fi 

l 
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CSEARCH A Any  command,  other  than  those  described  above,  transfers  the 

control  to  the  main  CLI  with  the  existing  offsets.  The  routine 
continues  the  spiral  search  by  the  CSEARCH  command  by  return- 
ing to  where  the  last  spiral  was  left  when  the  routine  was 
exited.  Any  other  commands  than  the  above  will  cause  an  EXIT 
of  this  routine,  but  with  the  last  specified  search  left  in 
progress.  This  command  allows  the  operator  to  re-enter  after 
leaving. 


2.5.2. 


Time  Changing  Commands 


Commands  under  this  subsection  allow  the  operator  to  change  the  time  for 
calculation  of  azimuth  and  elevation  angles  in  a number  of  ways  quite  similar 
to  the  actual  angle  changes  described  in  2.4.1. 


2. 5. 2.1.  MSLIP 


MSLIP  .jl  This  command  to  manually  change  time  is  analogous  to  MMOUNT, 

which  allows  manual  mount  control.  The  symbol  commands 
described  in  2. 4. 1.1.  under  MMOUNT  achieve  similar  effects  on 
time  offsets  under  MSLIP  as  they  did  on  angle  offsets  under 
MMOUNT.  (Not  valid  for  KALTRACK. ) Example: 
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MSLIP*^ 
SLIP  PER 


<0R  > ? (SEC.  ) The  system  asks  for  time  slip  in  seconds  per 

symbol  command'  ( < or>  ). 

Add  one  increment 
Subtract  one  increment 
Add  one 
Subtract  one 
Leave  this  routine 


TIME  OFFSET  IN  SEC  = X.XX 


2. 5. 2. 2.  CSLIP 


CSLIP  / 


This  command  for  continuous  change  in  time  is  analogous  to 
CMOUNT,  which  allows  continuous  mount  control.  The  symbol 
commands  described  in  2. 4. 1.2.  under  CMOUNT  achieve  similar 
effects  on  time  slip  rate  in  CSLIP  as  they  did  on  the  angle 
change  rate  under  CMOUNT.  Example: 


CSLIP^ 

SLIP  RATE?  (SEC/SEC) 

1 >* 

MAXIMUM  TIME  SLIP  (SEC) 

10; 

<  Jot  Oil-  Go  backward  in  time 

Speed  up  the  time  slip  rate  by  a factor  of  2 

< * Slow  down  the  time  slip  rate  by  a factor  of  2 

, * Go  one  third  as  fast 

/ Pause 

? Exit  and  give  total  offset  in  time 

TIME  OFFSET  IN  SEC  = X.XXXXX 


2. 5. 2. 3.  INCT 

INCT  y Increment  the  time  offset  by  a specified  amount  and  type  out 

the  total  offset.  (Not  valid  in  XALTRACK.)  Example: 

INCT^ 

INC  TIME  OFFSET  (H,M,S) 

0,4, -3  Operator  wants  a time  increment  of  4 minutes  and  -3  seconds. 

TIME  OFFSET  IN  HOURS  = .XXXXX  H 


r 


i 


2. 5. 2. 4.  SETT 

SETTil  New  specified  time  offset  is  set  and  the  previous  offset  is 

removed.  Example: 

SETT  J 

TIME  OFFSET  (H,M,S) 

1.2,3* 


2. 5. 2. 5. 

TIMEil 


2. 5. 2. 6. 

SDATE > 


SDATE/ 

ft  OF  DAYS  TO  GO  BACK  IN  TIME? 

2 / 

2. 5. 2.7.  SLIP 

SLIP/  This  and  subsequent  commands  may  be  considered  subsets  of 

CSLIP  since  what  they  achieve  can  be  done  equally  well  with 
CSLIP.  SLIP  initiates  a time  slip  as  in  CSLIP.  (Not  valid  in 
KALTRACK.)  Example: 

SLIP* 

SLIP  PATE?  (SEC/SEC) 

IP  Change  the  time  used  in  calculation  by  one  second  from  current- 

ly used  time  per  one  second  of  time. 

MAXIMUM  TIME  SLIP?  (SEC) 

10/ 


TIME 

This  command  simply  prints  out  the  current  time  from  the  clock 
without  any  offsets. 

SDATE 

This  command  allows  the  operator  to  go  back  in  time  by  typing 
a positive  response  to  the  system  query.  The  command  also 
clears  any  flags  that  have  been  mis-set  by  a time  glitch  from 
the  external  clock.  The  external  clock  often  gives  glitches 
which  cause  the  date  to  advance  by  one  day  per  glitch.  SDATE 
corrects  this  problem.  The  command  is  issued  after  the 
operator  has  determined  the  number  of  days  by  which  time  has 
been  shifted  due  to  glitches.  Such  a condition  generally 
causes  a runaway  mount  which  is  immediately  noticeable. 

Example : 
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2. 5. 2. 8.  PAUSE 


PAUSE  A 


2. 5. 2. 9. 

BACKjl 


BACKjl 

RATE 


hlOOjJ. 


2.6. 


This  command  pauses  a time  slip  already  in  progress.  (Not 
valid  in  KALTRACK.) 


BACKTRACK 


This  command  allows  a time  slip  in  progress  to  be  reversed  at 
operator-specified  speed  similar  to  BACK  in  the  spiral  search 
2. 5.1. 3.  (Not  valid  in  KALTRACK.)  Example: 


Go  back  at  the  present  rate 


ALPHABETICAL  INDEX  OF  COMMANDS 


» Vv.*'.  -•  • 
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No. 

Command 

Item 

1 

A 

2. 4. 1.1., 

2 

ABORT 

2. 5. 1.5. 

3 

BACK 

2. 5.1. 3. , 

4 

CARDINPUT 

CM 

CM 

CM 

CM 

5 

CMOUNT 

2. 4. 1.2. 

6 

CONT 

2. 2.2. 9. 

7 

CSEARCH 

2. 5. 1.6. 

8 

CSLIP 

2. 5. 2. 2. 

9 

E 

2. 4. 1.1., 

10 

EXIT 

2. 5. 1.4. 

11 

FILTER 

2. 2. 2. 4. 3 

12 

GETP 

2. 3. 2.1. 

13 

HEADER 

2. 3. 2.4. 

14 

INCM 

2. 4. 1.3. 

15 

INCS 

2. 4. 2.1. 

16 

INCT 

2. 5. 2. 3. 

17 

KALCARD 

2. 2. 2. 4.1 

18 

KALTRACK 

2. 2. 2. 4. 2 

19 

LIST 

2. 3. 2. 2. 

20 

MMOUNT 

2. 4. 1.1. 

21 

MSLIP 

2. 5. 2.1. 

22 

PAUSE 

2. 5.1. 2. , 

23 

PRINT 

2. 3. 2. 3. 

24 

QTRACK 

2. 2. 2. 5. 

25 

QUIT 

2. 2. 2. 6. 

26 

RDOS 

2. 2. 1.5. 

27 

READ 

2. 3. 1.2. 

28 

SATTRACK 

2. 2. 2. 4. 

29 

SDATE 

2. 5. 2. 6. 

30 

SEARCH 

2. 5.1.1. 

31 

SETM 

2. 4. 1.4. 

32 

SETS 

2. 4. 2. 2. 

2. 4. 1.2. , 2. 5. 2.1. , 2. 5. 2. 2. 


2. 4. 1.3. , 2. 5.2.1. , 2. 5. 2. 2. 


J 


33 

SETT 

34 

SLIP 

35 

START 

36 

STARTRACK 

37 

STOP 

38 

STORE 

39 

TIME 

40 

WRITE 

41 

XFER 

SYMBOL  COMMANDS 

42 

> 

2.4 

43 

< 

2.4 

44 

? 

2.4 

45 

2.4 

46 

9 

2.4 

47 

/ 

2.4 

48 

4- 

2.4 

49 

0 

2.4 

2. 5.2.4. 

2. 5. 2. 7. 

2. 2. 2. 7. 

2. 2. 2. 3. 

2. 2. 2. 8. 

2. 3. 1.4. 

2. 5. 2. 5. 
2. 3.1. 3. 
2. 3. 1.1. 


1.1.,  2. 4. 1.2.,  2. 5. 2.1.,  2. 5.2.2. 

1.1.,  2. 4. 1.2.,  2. 5. 2.1.,  2. 5. 2. 2. 

1.1. ,  2.4. 1.2.  , 2. 5. 2.1.  , 2. 5. 2. 2. 

1.1.,  2. 4. 1.2.,  2. 5. 2.1.,  2. 5. 2. 2. 

1.1.,  2. 4. 1.2.,  2. 5. 2.1.,  2. 5. 2. 2. 

1.2. , 2. 5. 2. 2. 

1.2. , 2. 5. 2. 2. 

1.2. ,  2. 5.2.2. 


SECTION  3 


PASSIVE  SATELLITE  TRACK 


3.1.  INTRODUCTION 


There  are  two  aspects  to  the  satellite  tracking  portion  of  the  coelostat 
project.  This  section  deals  with  the  passive  track  by  the  coelostat  mount. 

The  track  is  passive  in  the  sense  that  the  program  inputs  are  only  the  pre- 
specified set  of  orbital  parameters  for  an  epoch  and  no  real-time  information 
about  the  satellite  position  is  used  for  tracking.  The  sole  real-time  input 
for  passive  tracking  is  in  the  form  of  clock  time  and/or  angle  offsets  de- 
scribed in  the  User's  Manual.  These  offsets  serve  only  to  bring  the  satellite 
to  the  center  of  the  field  of  view  after  acquisition.  Active  satellite  track 
not  only  attempts  to  center  the  satellite  but  also  updates  the  orbital  ele- 
ments, and  is  described  in  Section  5. 

The  passive  satellite  tracking  program  utilizes  the  general  perturbation 
theory  first  proposed  by  Kozai  [4]  and  simplified  by  the  Space  Defense  Command 
(SDC).  The  Simplified  General  Perturbation  (SGP)  routine  calculates  the 
pointing  angles,  given  the  orbital  parameters  for  a specified  epoch  [3].  The 
perturbations  of  six  orbital  elements  of  a close  earth  satellite  moving  in  the 
gravitational  field  of  the  earth  without  air  resistance  are  derived  as  func- 
tions of  mean  orbital  elements  and  time.  No  assumptions  about  the  order  of 
magnitude  of  the  eccentricity  and  inclination  need  be  made,  i.e.,  the  trajec- 
tory solution  is  non-singular  even  at  zero  eccentricity  and/or  inclination. 
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Three  types  of  perturbations  have  to  be  taken  into  consideration:  1)  Secular, 

2)  Long-period  and  3)  Short-period. 

3.1.1.  Coordinate  Systems 

It  is  convenient  to  adopt  the  earth-centered  inertial  coordinate  system 
(ECI)  while  tracking  geocentric  satellites.  In  this  system,  the  earth's 
center  is  taken  to  be  the  origin.  The  x-axis  points  along  the  vernal  equinox, 
which  is  the  point  of  intersection  of  the  projection  of  apparent  path  of  sun 
on  the  earth's  surface,  and  the  equator  where  the  sun  goes  from  south  to  north 
in  its  apparent  annual  motion  along  the  ecliptic.  The  y-axis  lies  in  the 
equatorial  plane  and  is  perpendicular  to  the  x-axis  in  a direction  such  that 
the  z-axis,  which  completes  the  right-handed  system  points  towards  the  North 
Pole  (Figure  3-1).  Polar  coordinates  in  the  ECI  system  are  right  ascension 
and  declination  angles  as  shown  in  the  figure. 


To  define  the  orbital  coordinate  system  it  is  first  necessary  to  define 
the  orbital  parameters  of  a satellite  moving  in  an  elliptic  orbit  around  the 
earth. 


a semi-major  axis  of  the  elliptic  orbit  in  units  of  earth  radii  (ER). 
e eccentricity  of  the  orbit 

i inclination  of  the  orbital  plane  to  the  equatorial  plane 

ft  right  ascension  of  the  ascending  node,  i.e.,  the  point  where  the 


satellite  crosses  the  equatorial  plane  going  north  from  south. 


u 


argument  of  perigee,  i.e.,  the  angle  between  the  lines  joining  the 
earth  center  to  the  perigee  (line  of  apsides)  and  the  ascending 
node. 


E true  anomaly  is  the  angle  between  the  line  of  apsides  and  the  radius 

vector  pointing  to  the  current  satellite  position. 

In  the  orbital  coordinate  system,  axis-C  is  taken  along  the  line  of 
apsides  and  K -axis  is  perpendicular  to  the  orbital  plane  along  the  angular 
momentum  vector.  The  h-axis  is  also  in  the  orbital  plane  so  as  to  make  the 
system  right-handed  associated  with  the  unit  vectors  , i^  and  i^  (Figure 
3-2). 


Figure  3-2  Orbital  Coordinate  System 


The  transformation  of  coordinates  between  ^ ^ , £ and  the  x,y,z  axes  is 

affected  by  means  of  the  rotation  matrix  G . 

J o 


p __T  = G p = where 

K ECI  o H0C 


(3.1-1) 
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(3.1-2) 


G 

o 

C DC  to-  S S)S  oiCi 
S fiC  cj+  C S)S  wCi 

[ -C  nSuU>-  S ftC  coCi 
1 

1 -S  fts  QJ+  C £2C  coCi 
1 

S DSi 
-C  ftSi 

S u) Si 

i 

1 C uSi 

Ci  _ 

where  C0  = cos0  and  S0  = sin0  and  PECI  and  pQC  are  the  vectors  in  the  ECI  and 
Orbital  coordinates  respectively.  The  third  coordinate  system  which  will  be 
required  is  the  Local  Topocentric  System  at  the  observatory  site.  In  this 
system  the  E-axis  is  pointing  along  the  local  East,  the  N-axis  along  the  local 
North  and  the  Ze-axis  completes  the  right-handed  system  by  pointing  along  the 
local  zenith.  The  transformation  of  coordinates  between  the  x,  y,  z and  the 
E,  N,  Ze  axes  is  achieved  by  means  of  the  following  rotation  matrix  at 
sidereal  time  0 at  the  time  t, 


-S0 

C9 

0 

-S0C0 

-S0S0 

C0 

1 Px(t) 

(3.1-3) 

C0C0 

C0S0 

S0 

i.e.  , 0 = Greenwich  hour  angle  + local  latitude  (\  ) (Figure  3-3).  pT(t),  p^Ct) 
are  vectors  at  time  t in  the  topocentric  and  inertial  coordinates  respectively, 
and  0 is  the  latitude  of  the  site,  Figures  3-3,  3-4. 


3.2.  INITIALIZATION 


The  program  is  initialized  using  the  two-card  parameter  set  provided  by 
SPACETRACK  of  the  SDC.  The  format  of  data  received  is  described  in  Table  3-1 
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Table  3-1  (Continued) 
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=* — < 

and  contains  the  necessary  information  about  the  satellite  to  be  tracked, 
which  includes: 


1. 

epoch  day  and  fractional  day  of  the  year, 

2. 

angle  of  inclination  of  the  orbit,  i 

3. 

right  ascension  of  the  ascending  node,  ft 

4. 

argument  of  perigee,  to 

• 

5. 

mean  anomaly,  M 

6. 

mean  motion,  n,  in  revolutions  per  day 

7. 

first  time  derivative  of  the  mean  motion,  y;  and 

8. 

second  time  derivative  of  the  mean  motion,  jjr  . 

3.2.1.  Constants 

A number  of  constant  parameters  have  to  be  specified  either  as  data  or  in 
the  initialization  part  of  the  program. 

ag  Mean  equatorial  radius  of  earth  = 6378.145  km 

bg  Mean  polar  radius  of  earth  = 6356.7598  km 

f flattening  of  the  earth  = 3e  ^e  = 1 = 0.003352392 

a 298.25 

e 

GM  Geocentric  gravitational  constant  expressed  as  the  product  of 
the  universal  constant  of  gravitation  (G)  times  the  mass  of 
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20  3 -2 

the  earth  (M)  = 3.986012  x 10  cm  sec 
= 11467.87425  ER3DAY-2 

= 0.0055304177  ER3MIN~2 


3 -2  -2 

U Mass  function  of  earth  =1  ER  MIN 


k Geocentric  gravitational  constant  in  SPACEYRACK  format 

= 0.0743667785  ER3/2  MIN-1  = 0.001239446  ER3/2SEC_1 

3/2  -1 

= 107.0881622  ER  ' DAY 


Second  zonal  harmonic  coefficient  of  the  geopotential  function 
= 1082.28  + .3  x 10-6 


Third  zonal  harmonic  coefficient  of  the  geopotential  function 
= -2.562  x 10-6 


\ Geodetic  latitude  of  the  telescope  site  (Verona)  = 43.153291 


0 Longitude  of  the  telescope  site  (Verona)  = 75.625442W 

= -75.625442 


H Height  above  earth  mean  surface  at  the  telescope  (Verona) 
= 495.54'  or  2. 368095852E-5  ER 


. Vvl  ■.  i.  ■ • ‘ . *»•/;**  « - • 


The  card  data  supplies  the  mean  values  of  the  parameters  specified  in  the 
section  titled  Initialization.  These  parameters  remain  constant  till  they  are 
updated,  either  by  new  card-parameter  sets  or  by  active  tracking.  These  are 
the  mean  values  at  the  specified  epoch  time  and  are  subscripted  "0".  The  mean 
values  supplied  to  the  program  are  used  in  calculating  the  following  param- 
eters not  directly  specified  in  the  card  format. 


Lq  = mean  orbital  longitude 

= M + £2  + w (3.2-1) 

o o o 

The  mean  (Kozai)  semi-major  axis  is  calculated,  taking  into  account  the 
oblateness  of  the  earth.  For  a satellite  orbiting  around  a perfect  sphere 
with  an  isotropic  gravitational  field  we  have  the  following  relation  between 
the  semi-major  axis  of  the  ellipse  and  the  orbital  period  of  the  satellite: 

a3n2  = P (3.2-2) 

where  n is  the  mean  motion  in  revolutions  per  day  and  P=GM  is  the  gravita- 
. . . 3 -2 

tioi.al  constant  in  units  of  ER  DAY  . This  gives  the  semi-major  axis  ”a"  in 
units  of  earth  radii  too.  Due  to  the  oblateness  of  the  earth,  the  semi-major 
axis  given  by  (3.2-2)  has  certain  perturbation  introduced  into  it  which  is  a 
function  of  the  orbit  eccentricity  and  inclination. 


( 


j 





* 


a = a [1  + i <S-  i 62] 

O 0 0 


(3.2-3) 


where 


X A 3 

6 ' 2 


a2(l-e2)3/2 

o 


3.2.. 
(1  - -r  sin  l ) 

2 O 


(3.2-4) 


Now,  aQ  is  the  mean  value  of  the  semi-major  axis  in  units  of  earth  radii  at 
the  epoch,  corresponding  to  the  mean  epoch  values  supplied  by  the  card  param- 
eters. The  derivatives  of  ft  and  u have  been  expressed  as  functions  of  aQ, 
eQ  and  i so  that  they  result  in  imparting  continuously  varying  perturbations 
to  the  satellite  orbit,  and  are  constant  for  a given  set  of  parameters  at  an 
epoch. 


T ** 

2 cos  i (3.2-5) 

2,.  2.2  0 

ao(1_eo) 

2 T7  2TT  (2  - 2 Sln  1o)  (3.2-6) 

a(l-e) 
o o 

The  variation  of  ft  is  such  that  the  plane  of  the  orbit  rotates  about  the 
earth's  polar  axis  in  a direction  opposite  to  that  of  the  motion  of  the  sat- 
ellite with  the  mean  rate  of  rotation  given  by  Equation  3.2-5.  Similarly,  ci 
causes  the  line  of  apsides  to  rotate  about  the  satellites  angular  momentum 
vector.  It  can  be  checked  that  if  i >63°  26.7',  the  line  of  apsides  will 
regress,  while  if  i <63°  26.7',  it  will  advance  (Equation  3.2-6).  <L  is 
negative  in  the  former  case  and  positive  in  the  latter  case. 
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3.3. 


SECULAR  PERTURBATIONS 


As  explained  in  the  preceding  section,  and  impart  perturbations  to 

the  orbit  which  cause  rotations  of  the  orbital  plane  and  the  line  of  apsides 
about  the  polar  axis.  The  first-order  derivative  of  the  mean  motion,  n,  is 
mainly  the  effect  of  a decaying  orbit  and  is  also  responsible  for  secular 
perturbations  of  close  earth  satellites  which  experience  appreciable  drag 
effects.  Secular  perturbations  are  neither  functions  of  <*)  nor  M but  of  time. 
The  time  derivatives  J)0»  &Q  and  nQ  appearing  in  the  expressions  for  these 
perturbations  are  constants  for  a given  epoch  and  therefore  cause  perturba- 
tions to  be  monotonous  functions  of  time. 

At  - t - t = Time  since  epoch  (3.3-1) 

Since  the  mean  anomaly  M is  defined  as 

M = n(t  - t ) (3.3-2) 

o 

the  change  in  mean  anomaly  since  epoch  due  to  the  change  in  the  mean  motion  is: 
AM  = n At  + no  At2  + no  At3  (3.3-3) 

0 T “6 

Similarly,  J,  and  o cause  the  following  changes: 

o o 


Aw  = i^o  it  = change  in  argument  of  perigee  since  epoch;  (3.3-4) 
Aft  = ft^  At  = change  in  RA  of  ascending  node  since  epoch. 

The  intermediate  mean  values  of  the  parameters  which  will  be  required  for 
calculating  the  subsequent  long  and  short-period  perturbations  are  now  cal- 
culated. 

L = L + Am  + Aft  + Aw  (3.3-5) 

m o 

where  L is  the  mean  orbital  longitude  defined  in  Equation  3.2-1. 


The  equation  for  em  is  obtained  by  assuming  that  perigee  does  not  change  with 

time,  thus:  q = a (1-e  ) = a (1-e  ).  The  inclination  of  the  orbit  does  not 

^ o o m m 

suffer  any  secular  perturbations  [4]. 

3.4.  LONG  PERIODIC  PERTURBATIONS 


The  long  periodic  perturbations,  by  definition  [4],  are  functions  of  the 
argument  of  perigee.  Since  u is  a quantity  with  a period  of  variation  much 
larger  than  the  orbital  period  of  the  satellite,  the  perturbation  terms  which 
are  functions  of  go  are  said  to  be  long  periodic.  The  quantity  to  itself 
undergoes  secular  perturbations  due  to  ^ as  explained  in  the  preceding  sec- 
tion. Hence,  it  causes  the  line  of  apsides  to  completely  rotate  about  the 
satellite  angular  momentum  vector  in  a period  of  time  which  is  much  larger 
than  the  orbital  period.  Therefore,  perturbations  dependent  upon  to  are  said 
to  be  long  periodic. 


First,  e and  go  are  transformed  to  a new  set  of  variables 
m m 


A A 

a = e cos  to  : a = e sin  10 
xn  m m yn  mm 


(3.4-1) 


Then  the  effect  of  long  periodic  terms  is  calculated  as  follows: 


6 =A  tl 


a sin  i 
e m 

a / , 2 v 

m (1-e  ) 


(3.4-2) 
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mktri  *+ 


e cos  w 
m m 


e sin  w ^ . 

m m 2 L 


(3.4-3) 


- *2  > 1/2 

xnL  y"L 


(3.4-4) 


tan  yn. 


0 < u < 2 * (3.4-5) 

L 


The  restriction  of  to,  in  3.4-5  indicates  that  <o  must  be  calculated  as 

L It 

mod(  2 tt).  The  long  periodic  term  on  L gives 


, 1 jt  3 + 5 cos  i 

L - — o , a m 

m 4 L xn,  r : 

L 1 + cos  l 

m 


(3.4-6) 


0 < Ll  <2  TT 


Now,  Kepler's  equation  is  solved  for  the  eccentric  anomaly  using  the  function 
EXANM. 


E.  = EXANM  (L.-  to  - n , e ) which  solves  M = E - e sin  E by  simple 
L L L m L iteration 

(3.4-7) 


True  argument 
of  latitude: 


= VL+  WL 


V 

1/2 

H1_'l 

. 

tan  EL  + to 
2“  L 


(3.4-8a) 


iw  • * ;»  • a 


R.  = a ( l-er  cos  E.  ) 
L m L L 


( 3 . 4-8b ) 


Semi-latus  Rectum:  PT  = a (1-ef) 

Li  m L 


(3.4-8c) 


Transverse  com-  Rv  P 

ponent  of  velocity  — - — 


vector: 


( 3 . 4-8d ) 


Radial  component 

of  velocity  vector:  R =~\fu a"~  eL  sin  E 


( 3 . 4-8e ) 


where  v = true  anomaly. 


3.5.  SHORT  PERIODIC  PERTURBATIONS 


The  short  periodic  perturbations  are  functions  of  the  true  argument  of 
latitude  u,  which  itself  has  the  same  period  as  the  orbital  period;  hence,  the 
perturbations  are  said  to  be  short  periodic. 


, A 1 , e,2 

6S  - 4 J2  (p-} 

L 


( 3 . 5-la ) 


Range:  R = R.  + 6 (sin  i cos  2u)  PT 

^ L»  o m Li 


(3.5-lb: 


True  argument 

of  latitude:  uc  = u - ^ 6 (6-7  sin  i ) sin  2u 

a L i b m 


(3.5-lc) 


Inclination:  ic  = i +3  6 sin  i cos  i cos  2u 

s m S m m 


(3.5-ld) 


...  * / * - 


'■*  '*>  * 


1 


R.A.  of 

ascending  node:  fi„  = fi  + 3 6-  cos  i sin  2u 

b m b m 


(3.5-le) 


3.6.  OSCULATING  POSITION  AND  VELOCITY  VECTORS 

The  new  position  and  velocity  vectors  can  now  be  calculated  using  the 
orbital  parameters.  The  vector  along  the  radial  direction,  U_  is  given  as 
follows : 


U = u i + u i + uk 
x yJ  z 


(3.6-1) 


where 


u^  = cos  us  cos  fig  - sin  ug  sin  fig  cos  ig  (3.6-2a) 


u = cos  Ug  sin  fig  + sin  ug  cos  fig  cos  ig  (3.6-2b) 


u = sin  u„  sm 

Z b b 


( 3 . 6-2c ) 


The  unit  vector  V_,  perpendicular  to  the  radial  vector,  is  along  the  direction 
of  the  transverse  velocity  component  Rv  in  the  orbital  plane. 


V = -sin  u_  cos  fi_  - cos  u_  sin  ^ cos  i„ 
x S S S b S 


V = cos  u_  sin  i_ 

Z b b 
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( 3 . 6-3a ) 


V^  = -sin  Ug  sin  fig  + cos  Ug  cos  ftg  cos  ig  (3.6-3b) 


( 3 . 6-3c ) 


The  vector  W_  which  completes  the  triad  is  perpendicular  to  the  orbital  plane, 
i.e.,  parallel  to  the  direction  of  the  angular  momentum  vector. 


ux  = sin  is  sin 


(3.6-4a) 


= -sin  ig  cos  ns 


( 3 . 6-4b ) 


U)  = COS  1 
z S 


(3.6-4c) 


The  components  of  the  above  unit  vectors  are  along  the  earth-centered  inertial 
coordinate  system. 

Now,  the  radius  and  velocity  vectors  can  be  written  in  terms  of  these  vectors 
and  their  magnitudes. 


R = RgU 


( 3 . 6-5a ) 


R = R U + Rv  V 


(3.6  — 5b ) 


l5.|=  ^ t R2z)U2 


(3. 6-5c ) 


3.7.  CALCULATION  OF  POINTING  ANGLES 

Site  ,ector_:  The  site  vector  is  obtained  as  a function  of  the  sidereal  time 

The  rate  of  rotation  of  earth  is: 
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e 


0.00007292116  rad. /sec. 


The  local  sidereal  time  or  hour  angle  (Figure  3-3)  is 


®LHA 


®GHA 


Aw 


where  0 is  the  hour  angle  at  Greenwich  meridian  and  Aw  the  west  longitude 

brlA 

of  the  observer's  meridian  in  degrees.  If: 


0 = 99°. 6909833  + 36000°. 7689  T 

o u 


where 


_ _ J.D.  - 2415020.0  = tropical  centuries 

u = 36525 

since  noon  January  0,  1900  and  J.D.  = Julian  date  which  is  obtainable  from  a 
calendar,  then  the  Greenwich  hour  angle  is'  expressed  as 


6GHA(k)  = 


Atk 


d0 

dt 


From  this  the  recursion  equation  for  the  sidereal  times  at  k and  k+1  can  be 
derived  to  be: 


9(k+l) 


0,.  . + 7.292116  x 10-5  At, 

(k)  k 


where  At^  is  in  seconds,  and  At^  = - t^ 


m 


0 


Geodetic  latitude  of  the  site 


H 

f 

C 

S 


Height  above  sea  level 

Flatness  ratio  of  earth  = ae  be  = 1/298.25 

a 

e 


A 

A 


[1  - (2f-f2)  sin20]“1/2  ER 
C[l-( 2f-f 2 ) ] ER 


(3.7-1) 


where  ER  is  the  length  unit,  "Earth  radius".  Now,  the  site  vector  in  the 
earth-centered  inertial  coordinates  is  calculated. 


Rs(k+1) 


|(C+H)  cos0  cos9^  + 1^ 

! (C+H ) cos0  sin6.,  , , 

i (k+1), 

!(S+H)  sin0 


(3.7-2) 


The  vector  from  the  site  to  the  satellite  in  the  inertial  coordinates  is: 


o 


' (k+1) 


- Rs(k+1) 


(3.7-3) 


which  is  transformed  to  the  topocentric  system  by  the  transformation  matrix  G 
in  Equation  3.1-3. 


P 

E 

j 

-sin0 

i 

cose 

0 ~1 
1 

?T(k+l ) 

= 

% 

= 

-sin0cose 

-sin0sin6 

cos0 I t (k+1) 

RZe 

k+1 

cose  COS0 

cos0sin6  ^ 

| 

sin0 1 

— 1 k+1  (3.7-4) 

3-23 


r y 

The  components  P^.,  P^  and  PVg  point  towards  the  local  (site)  north,  east  and 
zenith  respectively. 

Calculation  of  azimuth  (A)  and  elevation  (h)  angles  from  the  vector  P^(k+1) 

I at  time  k+1  is  then  a trigonometric  procedure  (Figure  3-4). 

A = tan-1  ( P_/  PM)  (3.7-5) 

E N 

and 

II  I 

These  values  are  given  as  input  angles  to  the  mount  through  the  computer-mount 
interface. 

3.8.  DISCUSSIONS 

I 

The  passive  satellite  tracking  routine  was  used  successfully  to  track  a 
number  of  satellites.  Some  of  the  satellites  were  visible  in  the  viewing  tele- 
scope only  and  were  brought  to  the  center  of  the  field  of  view  using  the 
mount  control  routines.  Typically,  the  satellites  are  visible  for  no  more 
than  a few  minutes  (5-8  minutes  at  most)  and  sometimes  are  quite  faint.  A 
TV  monitor  must  have  a sensitive  camera  to  be  able  to  display  them.  Inflow 
of  its  position  measurements  to  the  active  tracking  program  must  also  be  as 
automated  as  possible  to  efficiently  utilize  the  short  time  available 
, for  observation. 


SECTION  4 


ACQUISITION 

4.1.  INTRODUCTION 

It  may  not  always  be  possible  to  obtain  up-to-date  values  of  satellite 
orbital  parameters,  or  the  error  in  the  parameters  may  be  just  large  enough  to 
cause  the  object  being  tracked  to  fall  outside  the  field  of  view  of  the  tele- 
scope. Therefore,  if  the  object  to  be  tracked  is  not  visible  in  the  scope,  it 
would  become  necessary  to  search  for  it.  Moreover,  the  nature  of  the  error  in 
angular  position  may  be  assumed  to  be  a near-constant  bias  (in  azimuth  and 
elevation,  see  Figure  4-1)  in  any  particular  pass  of  the  satellite.  Such 
error  appears  as  a constant  offset  from  the  direction  pointed  at  by  the  tele- 
scope, so  that  a search  can  be  conducted  while  the  telescope  is  in  the  track- 
ing mode.  A spiral  scan,  performed  while  the  "telescope  mean  direction" 
tracks  the  calculated  satellite  position,  is  the  search  method  selected. 

Spiral  search  has  some  obvious  advantages  as  compared  to  some  other  scan 
methods  used  in  radar  tracking;  the  methods  being  helical.  Palmer,  raster  or 
nodding  scan. 

The  spiral  radius  and  angular  position  are  both  taken  to  be  proportional 
to  "square  root  of  time  from  the  spiral  initialization."  The  "square-root-t" 
spiral,  as  it  will  be  called  henceforth,  was  chosen  because  it  attains  a 
constant  linear  velocity  as  the  angle  increases  as  opposed  to  the  spiral  where 
the  radius  and  angle  are  some  other  functions  of  time.  Since  the  radius  and 
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Figure  4-1  Nature  of  Tracking  Error 
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the  angle  are  both  taken  to  be  proportional  to  the  same  parameter,  namely 
"square-root-t" , an  Archimedes  spiral  is  obtained.  The  spiral  of  Archimedes 
(linear  relation  between  radius  and  angle)  gives  a uniformly  expanding  spiral 
and  all  the  desired  "area"  can  be  covered  at  an  operator-specified  speed.  A 
search  during  the  track  mode  will  consist  of  the  square-root-t  spiral  in  the 
LOS  (line-of-sight  coordinate  system  to  be  defined  in  the  next  section)  system 
and  will  appear  to  be  a stretched  and  distorted  spiral-like  motion  of  the  line 
of  sight. 

4.2.  SPIRAL  SEARCH  BY  THE  LINE  OF  SIGHT  (LOS) 


Searching  and  tracking  must  be  done  simultaneously,  which  is  to  say  that 
a spiral  search  about  the  calculated  satellite  position  is  to  be  performed 
while  the  mount  is  in  the  tracking  mode. 

Let  us  define  the  LOS  coordinates  such  that  the  y-axis  points  along  the 
LOS  (when  no  spiralling  is  being  attempted),  the  x-axis  is  perpendicular  to 
the  y-axis  and  parallel  to  the  elevation  mirror  rotation  axis.  The  z-axis 
completes  the  right-handed  triad.  When  no  search  is  being  attempted,  the  LOS 
lies  along  the  y-axis,  i.e.,  along  the  direction  specified  by  the  instanta- 
neous calculated  values  of  the  azimuth  and  elevation  angles.  The  plane  xz  is 
parallel  to  the  image  plane  with  z-axis  being  the  "vertical"  and  x-axis  the 
"horizontal"  (see  Figure  4-2). 

If  (u^,  u^_,  u^)  are  unit  vector  components  in  the  LOS  system,  then  the 
corresponding  vector  (u^.,  u^,  uzg)  in  the  local  topocentric  system  is  obtained 

* 
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Where  A and  h are  the  calculated  azimuth  and  elevation  angles,  respectively. 


T 

Let  UL0S  ” (V  Uy ’ uz^  be  defined  as  a unit  vector  along  the  true  LOS 
and  Uxz  its  Pr°3ection  on  the  xz  plane.  The  tip  of  the  vector  uxz  is  required 
to  describe  an  operator-specified  spiral  in  the  xz  plane.  The  magnitude  of 
uxz,  p = ' uxz  1 must  be  allowed  to  grow  as  the  arg(uxz>=0  increases. 


We  have 


ux  = pcos0  ; uz  = Psine 

uy  = n'I'VuJ  = Jl- p2  (4.2-2) 


where  p and  0 are  calculated  as  functions  of  time  according  to  the  specified 
square-root-t  spiral.  After  transforming  uLQS  to  the  topocentric  system  uT> 
the  new  azimuth  and  elevation  angles  are  computed  from  the  following  expres- 
sions . 


The  angles  (A',  h')  are  supplied  to  the  mount,  through  the  interface,  as  the 
pointing  angles  during  the  track-search  mode,  instead  of  the  original  A and  h 
angles. 


4.3.  SPIRAL  GENERATION 

4.3.1.  Spiral  Parameters 

A spiral  is  specified  by  the  relation  P = f ( © > or  by  specifying  p and  9 
as  functions  of  some  common  independent  parameter.  Now,  spiral  time  is  de- 
fined as  "the  time  since  the  initiation  of  the  spiral." 

Spiral  time  = t = Current  time  -(Time  at  p = 0 = o)  (4.3-1) 

p = k^ft  ; 9 = k2>Tt  (4.3-2) 

p = ^1_  9 : Archimedes  Spiral  (4.3-3) 

k2 

dp  = kl  d£  = k2  = ^2_  d£_  (4.3-4) 

dt  2 rt  ’ dt  2jt  k^  dt 

Since  the  spiral  is  initiated  with  the  object  of  searching,  it  will  be  ad- 
vantageous to  maintain  constant  the  rate  at  which  the  center  of  the  field  of 
view  moves  along  the  spiral.  This  statement  and  the  advantages  of  a square- 
root-t  spiral  will  become  clear  presently. 
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The  length  of  an  element  of  the  spiral  between  ( p , 6)  and  (P  + d P , 6 


From  Equation  4.3-6  it  is  clear  that  the  rate  of  change  of  the  spiral  length  s 
approaches  a constant  with  increase  in  time.  If  P and  6 were  proportional  to 

r — dS  • 

t instead  ofjt,  then  would  have  increased  rapidly,  resulting  in  a decreas- 
ing "view  time"  and  making  viewing  and  object  identification  impractical  after 
a certain  time. 


4.3.2. 


4. 3. 2.1.  View  Time 


The  view  time  is  defined  as  the  time  for  which  an  object  passing  through 

the  aperture  center  would  remain  in  the  field  of  view,  and  is  determined  by 

ds  • • • • 

the  spiral  speed,  The  actual  length  of  the  field  of  view  is  a complicated 


function  of  the  aperture  diameter  (AD)  and  the  spiral  parameters.  However, 


for  0 > 2 it  , the  described  spiral  can  be  assumed  to  be  fairly  straight  in 
parts , so  that 


Length  of  field  of  view  = AD  , and 

Time  of  view  = T : Operator-specified  quantity. 


ds  • 

Let  T be  the  view  time  at  0 = 4ir  , when  has  become  fairly  constant. 


ds 


vdt/ 


0=47T 


klk2  I 1 


16ir 


+ 1 = 


~ AD 
T 


(4.3-7) 


klk2. 


16tt  +1 


8tt 


4. 3. 2. 2.  Overlap  Fraction 


The  overlap  fraction  specified  by  the  operator  refers  to  the  overlap  of 
the  apertures  at  9 and  0+2tt  . The  quantity  is  to  be  operator-specified  and  is 
the  increase  in  the  radius  p as  0 increases  by  2tt  , in  units  of  aperture 
diameter.  From  Equation  4.3-3  we  have 

P0+2tt  P0  = 2*  = f.AD  (4.3-8) 

k2 

where  f is  the  operator-specified  fraction.  Equations  4.3-7  and  4.3-8  can  now 
be  solved  for  the  two  unknowns  k^  and  k^.  The  aperture  diameter  is  obtained 
from  the  angular  field  of  view  of  the  telescope,  which  is  equal  to  5 arc 
minutes.  The  AD  for  the  unit  vector  considered  in  section  4.2.  will  be 
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AD 


angle  x length  from  origin  along  LOS 


Lj 


, 5 TT 

( 60  X 180  ’ 1 


(4.3-9) 


Using  Equations  4.3-7  through  4.3-9  gives: 


-4  -1  1/2 

k = 8.1928652  x 10  (fT  ) 


k2  = 3.5393177  (Tf) 


-1/2 


(4.3-10) 


The  third  quantity  that  the  operator  must  specify  is  the  maximum  angular 


deviation,  p , to  be  searched  in  arc  minutes.  The  program  will  calculate 
’ max 


and  output  the  search  time. 


_ P 


max 


max 


, P 

( max 


60 


t ,2  1 

180°  72  S6C- 

kl 


(4.3-11) 


4.4. 


OPTIONS  AND  CONTROL  COMMANDS 


4.4.1.  Options 


A number  of  options  are  included  to  aid  in  a search.  The  initialization 
options  allow  the  operator  to  specify  the  following  quantities. 


a.  The  view  time,  which  is  the  time  taken  by  the  aperture  to  traverse 
one  aperture  diameter  along  the  spiral,  is  also  the  approximate  time 


4-10 


„ , VvT-.  •$.  — 4fVl  k-  • • * - 


J 


r 


that  a point  will  stay  in  the  field  of  view.  A small  view  time  will  ‘ 

result  in  a fast  spiral  and  vice  versa.  ' 

b.  The  overlap  fraction  determines  the  compactness  of  the  spiral  and 
should  be  less  than  or  equal  to  unity  if  no  area  around  the  cal- 
culated point  is  to  be  missed.  Too  small  a value  of  f will  result 

in  a larger  search  time  t 

° max 

c.  Maximum  search  angle  is  determined  by  the  expected  error  offset,  and 
a value  of  2-6  times  the  angular  field  of  view  should  be  quite 
suitable. 

d.  The  search  stops  when  Pmax  is  reached  and  the  mount  continues  along 
the  latest  offsets  till  a command  to  either  backtrack  or  abort  the 
search  is  issued. 

e.  The  backtrack  command  retraces  the  spiral  backwards  at  an  operator- 

specified  speed.  It  may  be  issued  after  the  p is  reached  or 

max 

after  a pause  command.  The  spiral  can  be  stopped  as  soon  as  the 
object  is  viewed  and  the  mount  will  continue  at  the  latest  offsets. 

If  the  aperture  overshot  the  object,  backtrack  will  be  necessary. 

The  operator  may  induce  a slower  backtrack  so  as  not  to  miss  the 
object  again,  and  stop  it  as  soon  as  the  object  reappears  in  the 
field  of  view. 


f 


4.4.2. 


Commands 


The  satellite  track  program  has  two  command  line  interpreters  (CLI).  The 
main  CLI  deals  with  the  total  tracking  program  and  the  search  routine  CLI 
operates  only  when  the  program  is  in  a search  mode.  At  least  the  first  four 
letters  of  the  commands  must  be  entered. 

SEARCH : This  command  transfers  the  control  to  the  search  CLI  even  while  the 

main  program  continues  to  track.  The  computer  responds  by  asking  for  values 
of  view  time,  aperture  percentage  overlay  desired  and  the  maximum  angle  of 
search.  The  operator  is  expected  to  input  the  view  time  in  seconds,  the 
fraction  of  the  total  aperture  not  overlapping  after  the  spiral  goes  through  2ir 
in  terms  of  percentage  of  total  aperture,  and  the  maximum  angle  of  search  in 
arc  minutes.  When  all  the  three  questions  are  answered  correctly  by  the 
operator,  the  program  starts  the  spiral  after  putting  out  the  expected  time  of 
search,  depending  on  the  maximum  angle  of  search. 

PAUSE : When  the  object  is  sighted,  or  for  any  other  reason,  the  operator  can 

issue  this  command  to  interrupt  the  search.  The  track  continues  by  adding  the 
latest  offsets  to  the  calculated  azimuth  and  elevation  values. 

BACK:  If  the  object  is  sighted  during  the  search  and  it  overshoots  the  center 

of  field  of  view  before  issuing  the  PAUSE  command,  BACK  command  can  be  used  to 
backtrack  along  the  same  spiral.  The  computer  asks  for  the  rate  of  backtrack- 
ing in  terms  of  percentage  of  the  rate  at  which  the  original  spiral  was  being 


described.  An  operator  response  of  a value  between  0 and  ICO  causes  a 
slower  backtrack  and  a negative  value  causes  the  spiral  to  continue  in 
the  forward  direction  again. 

EXIT : The  EXIT  command  takes  the  task  out  of  search-routine  with  the  current 

azimuth  and  elevation  offsets. 

ABORT : This  command  aborts  the  search  and  zeros  the  existing  azimuth  and 

elevation  offsets  introduced  by  the  search. 

: Any  other  command  transfers  the  control  to  the  main  CLI . 


CSEARCH : This  command  brings  the  control  back  to  the  search  CLI. 

4.5.  RESULTS 


The  search  routine  has  proved  to  be  a valuable  tool  during  tracking 
experiments.  It  has  also  been  used  to  determine  the  fields  of  view  of  the 
main  telescope  and  the  viewing  telescope  on  the  television  screen.  By  deter- 
mining this  correspondence,  it  is  possible  to  observe  the  object  on  the  TV 
screen  and  say  whether  it  will  be  visible  in  any  or  both  of  the  telescopes  or 
not.  It  has  also  been  used  to  determine  the  portion  of  the  photographic  film 
that  is  exposed  through  the  main  telescope.  It  has  rarely  been  necessary  to 
call  the  search  routine  to  serve  the  purpose  for  which  it  has  been  designed. 
Stars  are  always  visible  in  the  center  of  the  field  of  view  because  their  RA 
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and  DEC  values  are  accurately  tabulated  in  [3].  Time  duration  for  which  a 
satellite  is  above  the  horizon  is  small;  however,  most  satellites  were  visible 
in  the  viewing  telescope  at  the  first  attempt.  Also,  some  of  the  satellites 
are  either  extremely  dim  and  therefore  barely  visible,  or  blinking  periodi- 
cally due  to  spinning.  Hence,  it  will  be  difficult  for  the  viewer  to  make  a 
decision  to  search  for  a satellite  if  it  is  not  visible  in  the  viewing  scope. 
The  invisibility  could  be  attributed  to  large  blinking  period,  small  amount  of 
reflected  light  or  erroneous  orbital  parameters.  Only  in  the  latter  case 
would  it  be  advisable  to  search,  and  even  then,  sighting  an  acquired  satellite 
cannot  be  guaranteed. 


Any  misalignment  of  the  mount  can  be  detected  using  the  search.  If  a 
star  is  not  visible  in  the  center  of  the  scope,  the  search  routine  can  be  used 
to  acquire  it.  The  offsets  required  to  do  so  will  be  the  mount  misalignments. 
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SECTION  5 

ACTIVE  SATELLITE  TRACK 

5.1.  INTRODUCTION 

The  passive  satellite  tracking  routine  is  fed  orbital  elements  of  the 
satellite  from  the  card  parameters  which  have  been  calculated  for  some  speci- 
fied epoch.  These  parameters  could  have  been  calculated  a few  days  earlier 
and  therefore  may  differ  slightly  or  substantially  from  the  current  values. 
After  the  satellite  is  acquired  in  the  field  of  view,  a number  of  measurements 
of  its  true  azimuth  and  elevation  can  be  obtained.  These  measurements  can  be 
used  either  on-line  or  off-line  to  update  the  orbital  parameters. 

The  active  tracking  program  was  designed  to  take  maximum  advantage  from 
the  existing  orbit  prediction  SGP  routine  described  in  the  section  on  passive 
tracking.  A model  of  the  orbital  parameter  variation,  based  on  the  secular 
perturbations,  was  formulated.  The  filter  is  designed  sc  as  to  accept  azimuth 
and  elevation  angles  as  the  observations. 

5.2.  LINEARIZED  DISCRETE  KALMAN  FILTER 


The  theory  of  Kalman's  recursive  minimum  variance  filter  has  been  des- 
cribed in  innumerable  papers[5]  and  texts[6].  Only  a summary  of  steps 
required  will  be  presented  here. 
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There  is  a nonlinear  relation  between  the  states  at  adjacent  time  in- 
stants and  also  between  the  observations  and  the  states  (which  in  this  case 
will  be  the  orbital  parameters). 


System  dynamics 


X.  . = f(X,  ) + 

k+1  k k 


(5.2-1) 


Observation  equation 


= MX. 

k+1  k+1  k 


(5.2-2) 


Where  X,  is  the  six-dimensional  state  vector  and  Z,  , is  the  two-dimensional 

k k+1 

observation  vector.  The  nonlinear  functions  f(X^)  and  h(X^)  are  six-dimen- 
sional and  two-dimensional  vectors,  respectively.  In  case  of  linear  systems, 
they  reduce  to  and  H^X^,  where  and  are  (6x6)  and  (2x6)  matrices, 

respectively  and  independent  of  the  state  vector  X^. 


The  vector  is  the  system  noise  representing  the  uncertainty  in  the 
system  model  and  is  assumed  to  be  Gaussian-distributed  with  zero  mean  and 
covariance  Q^.  The  measurement  noise  is  also  assumed  to  be  Gaussian- 


distributed  with  zero  mean  and  covariance 


V 


E(Wk)  = [o]6  ; E(Uk)  = [o]2 


E((W  = Vk*  ; E(W  = Vk* 


(5.2-3) 


Matrices  Qk  and  must  be  positive  semidefinite  and  definite,  respectively. 
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Here,  6 ^ is  the  Kroenecker  delta  function,  emphasizing  that  there  is  no 

correlation  in  time  in  the  noise  sequences. 

To  initiate  the  filter,  an  initial  estimate  of  the  state,  X and  its 

’ o 

covariance,  is  required.  This  may  be  obtained  from  the  preliminary  know- 
ledge about  the  orbit. 


The  predicted  value  of  the  state  at  a stage  k+1  is  obtained  from  the  best 
prior  estimate  at  stage  k,  X^, 

Xk,l  = f(xk) 


and  the  covariance  prediction  is  obtained  using  the  linearized  state  equation, 
which  gives 


k 


if.  . 


i,  j = 1,. • .6 


k+1 


(5.2-5) 


where  the  (6x6)  matrix  is  a linearized  state  matrix.  If  the  predicted 
covariance  is  denoted  by  M^+^,  then 


M = E 
k+1 


{(X. 


k+1 


xk,i>  "W  - W } s 


F.P,  F‘ 
k k k 


+ Q,  (5.2-6) 


The  suboptimal  Kalman  gain  for  the  linearized  system  is: 


K.  , - M,  .H*  (H,  ,M.  ,H.  , 

K+1  k+1  k+1  k+1  k+1  k+1 


(5.2-7) 
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where  M,  , is  the  (2x6)  linearized  observation  matrix. 
k+1 


3h, 


k+1 


~3x 


= Ch. . ] 


k+i  / 7 

Tc+l 


i = 1,2 

ij'y  j = 1....6 

Xk+1 


(5.2-8) 


Both  F,  and  H,  , have  been  linearized  about  the  most  recent  estimate  of  the 
k k+1 

state,  which  in  this  case  is  X^+^.  A residual  is  formed  by  differencing  the 
observation  and  the  predicted  observation. 


yk+i = zk+i  - h(W  = residual  (5-2-9) 

The  final  corrected  estimate  is: 

Vi  = xk.i  * ww  (5-2-10) 

which  has  an  error  covariance  matrix  given  by: 

pk.i  ■ (I  - lWWmk.i  <5-!-n) 

where  I is  the  (6x6)  unity  matrix.  An  equivalent  and  numerically  more  attrac- 
tive form  of  Equation  (5.2-11)  is: 

Vi  * ” - Kk.i"k.i>“k.i(I  - "kAi’1  * ‘WW&i  <5-2-i2> 

5.3.  SYSTEM  DYNAMICS 


Six  quantities  are  required  to  specify  the  position  of  a satellite  in  the 
ECI  coordinate  system.  The  six  orbital  parameters  (or  some  of  their 
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combinations)  or  the  radius  and  radial  velocity  vectors  are  two  of  the 
common  choices.  The  state  vector  to  be  estimated  was  chosen  to  be  the  set 


of  six  parameters  because  they  undergo  perturbation  and  have  to  be  corrected 
as  the  observations  are  being  taken.  The  state  vector  is  therefore: 


X^  = (a  e i ft  w E)^ 


(5.3-1) 


where  the  subscript  k refers  to  time  t^.  The  variation  of  this  vector  as  a 
function  of  time  is  the  system  dynamics  equation  which  is  to  be  specified. 


5.3.1.  Recursion  Equation 


Since  an  = p,  and  n is  changing  with  time  by  virtue  of  n,  we  can 
write  the  time  derivative  of  a.  This  relation  will  be  used  often  in  this 
section,  whenever  a derivative  of  a function  with  respect  to  'a'  is  required 
and  the  function  contains  the  term  'n'  - the  mean  motion. 


da  . 1/3, _2.  -5/3 - 

37  ■ 11  (~)n  n 


ak+l  = ak  + akA  tk 


4 k n . 

ak  ‘ 3 n,  2 A k 

k 


a (1  . 4 £.  _AJ_) 

kU  3 2 n,  } 

k 


(5.3-2) 


r 


Assuming  that  the  perigee 


distance  does  not  change  with  time,  we  have 


q = a - ae  = a(l  - e)  = Constant. 


(5.3-3) 


Taking  the  derivative  of  Equation  (5.3-3) 


de  _ 1 - e da 
dt  = a dt 


l-e. 


ek  = ak 


*k 


(l-e,) 


4_  n ~k 
3 2 n. 


• l-e. 

, 4_  n k 

* ek+l  = ek  ' ’ 


(— ) - 
^ 3 1 2 


3 2 n, 


At, 


(5.3-4) 


Even  though  the  inclination  has  short  period  perturbations,  its  mean  value 
remains  constant. 


• Vl  = Xk 


(5.3-5) 


The  time  derivatives  of  ft  and  4)  have  been  discussed  in  the  Section  3-2-2. 
of  passive  tracking  and  they  are: 


COS  1, 
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r 


= + ir  J 


5 2.  , 

2 “2  _2, , _2 n2  U ' 2 n Xk' 


VHP 


“ k+l  = k + ^ k Atk 


(5.3-6) 


and  u) 


k+1  = ^ 


+ a)  At, 


(5.3-7) 


The  recursive  equation  for  the  eccentric  anomaly  can  be  derived  from  Kepler's 
equat ion. 


M = E - esinE  = n(t-t  ) 
o 


(5.3-8) 


Ek+1  - Ek  - e(sinEktl  ' sinEk)  = nk  Atk 


AEk  ” e[sinEk  + AE^cosE^  - sinE^]  = nkA  tk  (5.3-9) 


The  above  expansion  is  obtained  under  the  assumption  that  A t,  , and  therefore. 


AEkj  are  small. 


AEk  - ecos  Ek  . AEk  = nk  Atk 


Ek+1  Ek  + (1-ecos  E,  ) Atk 


(5.3-10) 


Equations  (5.3-2  through  5.3-7  and  5.3-10)  give  the  recursion  of  the  state 
defined  by  Equation  (5.3-1).  These  equations  are  used  to  get  the  predicted 
state  Xk+^  from  the  estimated  state  Xk  as  in  Equation  (5.2-u). 
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5.3.2.  Linearized  System  Matrix 

Since  there  is  a nonlinear  relation  between  X^+1  and  X^,  linear  Kalman 
filtering  techniques  cannot  be  directly  applied.  An  approximation  must  be 
made  and  the  first-order  sensitivity  matrix  of  X^+1  with  respect  to  X^  must 


be  derived[6] 


f 'I 

9f(V 


where 


(5.3-12) 


From  Equation  (5.3-1),  we  have 


Xk  = ak  ; Xk  = ek  ’ Xk  \ 1 


\ ~ ‘ Xk  wk  and  Xk  = Ek 


(5.3-13) 


1st  Row:  Semi -major  axis 


3fk  u ; Atk 

11  - 3ak  - 3 2 nk 


f,  . = 0 for  j = 2 - 6 

i: 


(5.3-14) 


r 


2nd  Row:  Eccentricity 


; {1-eJ 

n k . 

nk  *k  k 


• At 

= i + J 

3 2 n, 

k 


f23  = f24  = f 25  = f26  " 0 


(5.3-15) 


3rd  Row:  Inclination 


Since  mean  inclination  is  constant,  all  the  terms  in  the  third 
matrix,  except  the  (3x3)  term,  are  zero. 


row  of  the 


f3 

f = JJi  = 1 
33  3 \ 1 


f31  " f32  " f34  = f 35  = f36  = 0 


(5.3-16) 


-~'h  Row'  Right  ascension  of  the  ascending  node  is  a function  of  all  param- 
eters except  for  the  argument  of  the  perigee  and  the  eccentric  anomaly. 


f45  fU6  = 0 5 f44  = 1 


"41  = = a/Cl-e?)2  C°S  **  ^ 

K k 


C°s  ik  Atk 


2 2 , 2,3 

ak(1"ek) 


= + I J2  2~  k '2~2  Sin  *k  Atk 

ak(1-ek} 


(5.3-17) 


5th  Row:  Argument  of  Perigee  is  a function  of  all  parameters  except  ft  k and 


f,,  = fcc  = 0.  Also  fcc.  = 1 
54  56  55 


3fk  21  T nk  5 . 2 . , 

r51  3a,  " 4 J2  3 , 2,2  (2  " 2 Sin  1k)  A k 

k ak  (l-ek) 


3fk  _ T nk  ®k  5 2 . , 

f 52  3e,  ' 6 J2  2/  2,3  (2  " 2 Sln  1k)  Atk 

k ak<l-ek> 


f53  = ” J2  2~  2,2  Sin  Ak  COS  *k 

*k  (1"ek) 


(5.3-18) 


6th  Row:  The  eccentric  anomaly  is  a function  of  ak  and  ek  only,  so  that 


f63  = f64  = f65  = °' 


3 nk  Atk 


2 ak  (i“ekcos 


62  3e, 


n,  cos  E, 

= + (l-ek  cos  Ek)2  Atk 


(5.3-19) 


3f,  n ek  sin  E 

f = =1 — _ — At 

66  3E,  1 P ,2  \ 

k (1-e^cos  E^) 


% Vs.- . • A 
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Equations  (5.3-14)  through  (5.3-19)  define  the  (6x6)  system  dynamics 

matrix  which  was  introduced  in  Equation  (5.3-11).  The  system  noise 

matrix  defined  in  Equation  (5.2-3)  is  determined  experimentally  by  trying 

different  values  in  simulated  missions.  The  system  matrix  F^  is  calculated 

at  Xk+1;  i.e.,  the  predicted  components  of  Xk+1  = (i^,  Tk+1>  ^k+1, 

u)  , , , E,  , )^  are  used  to  calculate  the  various  elements  of  the  system 

k+1*  k+1  J 

matrix.  The  calculation  of  Fk  is  simplified  by  the  fact  that  out  of  a 
possible  36  terms,  only  twelve  are  non- zero  terms  which  require  recalcula- 
tion at  every  computation  stage,  and  three  more  are  equal  to  unity. 

The  state  matrix  F,  is  used  to  calculate  the  predicted  covariance  M,  , 
k k+1 

according  to  Equation  (5.2-6)  from  the  error  covariance  P,  at  the  preceding 

K 

stage. 

5.4.  OBSERVATIONS 


The  azimuth  and  elevation  angles  will  be  the  observations  for  the 
updating  of  the  orbital  elements.  There  is  a non-linear  relationship  between 
the  observations  and  the  orbital  elements  of  the  form  given  by  Equation 
(5.2-2)  which  must  be  explicitly  determined  to  be  able  to  derive  the  line- 
arized observation  matrix  given  by  Equation  (5.2-8). 

5.4.1.  Relation  Between  Angle  Measurements  and  the  Orbital  Elements 

The  radius  vector  to  the  satellite  in  the  orbital  coordinates  is  given 


by: 


(5.4-D 


' k+1 


ak+l(cOS  Ek  + 1 -ek+l) 


/ 1 2 vl/2  . _ 

ak+l(1  - ek+l)  Sln  Ek+1 


k+1 


This  vector  is  now  transformed  to  the  ECI  coordinates  using  the  transformation 


matrix  Gq  from  Equation  (3.1-1) 


G = Go  G.C 
0 1 


where 


[cos  Q 


k+1 


sin  ft 
0 


G. 

l 


k+1 


‘s 


(5.4-2) 


-sinft 

cosft 

0 


i 

i 

i 

> k+1 


( 5.4-3a) 


0 

cos  i 
sin  i 


0 

-sin  i 
cos  i 


( 5. 4-3b ) 


k+1 
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along  the  local  east,  north,  and  zenith.  The  measured  angles  are  then: 


The  actual  measured  angles  would  be  the  true  angles  given  by  this  equation, 
contaminated  by  the  instrument  measurement  errors  which  are  represented  by 
Equation  (5.2-2). 

5.4.2.  Linearized  Observation  Matrix 


In  Equation  (5.4-5),  p^k+l)  is  the  site-satellite  vector  in  the  ECI 
coordinates.  The  derivative  of  the  expression  for  h(X^  ) with  respect  to 
has  to  be  taken  by  the  application  of  the  chain  rule. 


3h(Xk+i) 


3X. 


k+1 


3h 

3Pj(k+l) 


3?I(k+1)  _ 3h 

3Xk+1  ' 3<^(k+l) 


3_ 

3 X 


• ^'Vk+l 

k+1 


(5.4-7) 


Considering  that  R<,  is  independent  of  X^+1,  since  it  is  the  site  vector. 
Equation  (3.7-2),  and  also  suppressing  the  subscript  (k+1)  for  simplicity, 
we  get: 


3h(X)  3h 
3X  " 3Pj 


3R  s 3h  8pT 

3X  ” 3ot  ' 3pj 
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= H1H2 


(5.4-8) 


where , 


A 9h 

H1  = -JL' 
1 3p. 


a (2x3)  matrix 


. A 9r 
^2  = "3x 


= a (3x6)  matrix 


so  that  H is  a (2x6)  matrix. 


5. 4. 2.1.  Calculation  of  = [hl^] 


The  derivative  of  the  vector  h,  with  respect  to  the  vector  Pj,  is 
obtained  simply  by  considering  that  the  derivative  of  arc  tangent: 


d . -1  . 1 du 

-7—  (tan  u)  = ^ 

dx  , 2 dx 

1 + u 


First  row: 


d°X 


(5.4-9) 


where  p , p , p are  the  components  of  p in  the  EC  I system.  The  ex- 

x y z j J 

pressions  for  ( p_  pXT  p ) in  terms  of  ( p P P ) are  obtained  from 
L N Ze  x y z 

Equation  (3.7-4). 
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(5.4-lCa) 


Proceeding  in  a similar  fashion,  other  terms  of  the  matrix  H'  can  be  derived 
to  give  the  following  expressions. 
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(5.4-lOc) 


Second  row: 
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Similarly: 
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EN 
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(5.4-lCf) 
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Equations  (5.4-10)  define  the  matrix  completely. 

5. 4. 2. 2.  Calculation  of  H,, 

The  Jacobian  of  the  transformation  (5.4-4)  is  the  matrix  which  is 
obtained  by  writing  explicitly  in  terms  of  the  orbital  parameters 

[Equations  (5.4-1)  through  (5.4-4)]  and  then  taking  appropriate  partial 
derivatives.  Since  the  relation  is  available  as  a matrix-vector  product  in 
Equation  (5.4-4),  can  be  calculated  columnwise  by  taking  derivative  of 
with  respect  to  each  of  the  orbital  elements. 


1st  Column:  Only  rg  is  a function  of  the  semi-major  axis  a. 


3r 

-4—  = G -r  = G 

da  o 3a  o 


cos  E - e 


i£  2 


e sin  E 


(5.4-lla) 


2nd  Column:  Only  rg  is  a function  of  the  eccentricity  e. 
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3rd  Column: 
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4th  Column: 
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5th  Column: 
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6th  Column: 
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Now,  the  linearized  measurement  matrix  H is  given  by  postmultiplying  by 
Equation  (5.4-8). 


The  measurement  error  covariance  matrix  R is  taken  to  be  positive 
definite  and  compatible  with  the  expected  error  variance  in  the  angle- 
measuring instruments.  Generally,  it  is  taken  to  be  a diagonal  matrix. 

Then  the  Kalman  gain,  the  corrected  estimate  of  the  six  orbital  parameters 
and  the  updated  covariance  matrix  are  computed  according  to  Equations  (5.2- 
7),  (5.2-10)  and  (5.2-11)  respectively.  The  whole  process  described  in 
Sections  5.3  and  5.4  is  repeated  at  each  observation  sampling  stage. 
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SECTION  6 

PROGRAM  DESCRIPTION 


6.1.  INTRODUCTION 

This  section  of  the  report  describes  the  computer  program  and  hardware 
interface  developed  by  PAR  for  driving  the  optical  tracking  mount  at  the  PATS 
tower  of  the  Verona  test  annex  of  the  RADC.  The  theory  supporting  the  track- 
ing programs  has  been  described  in  Sections  3,  4,  and  5 of  this  report.  The 
current  section  gives  a brief  description  of  the  FORTRAN  coded  programs,  at 
varying  levels  of  detail,  to  supplement  the  elaborate  documentation  provided 
in  the  computer  source  listing  in  the  form  of  comments.  The  listing  of  each 
routine  includes  a definition  of  each  of  the  internal  and  external  variables, 
common  blocks,  required  inputs  and  available  outputs,  detailed  description  of 
the  subroutine  function,  operations  performed  in  different  sections  of  the 
program  and  external  subprograms. 

The  system  has  been  described  with  supporting  flowcharts  of  varying 
complexities  from  overall  system  blocks  to  individual  subroutine  descriptions. 
Then,  a short  description  of  the  hardware  is  followed  by  an  index  of  programs. 

The  main  functions  that  the  overall  system  can  perform  are: 

a.  Track  stars  given  the  right  ascension  and  declination  from  star 
catalogs ; 

b.  Track  satellites  given  NORAD  satellite  card  parameters  using  a 
passive  tracking  routine; 

c.  Track  a satellite  and  update  its  parameters  by  filtering  its 
observed  azimuth  and  elevation  positions  using  the  Kalman  filtering 
technique . 

All  tracking  systems  have  a hard-wired  clock  pulse  input  for  the  Green- 
wich Mean  time  and  other  inputs  required,  depending  upon  the  tracking  program 
in  control.  The  STARTRACK  program  requires  the  right  ascension,  declination, 
brightness  magnitude  and  a catalog  number  for  the  star  to  be  tracked  for  its 
STARLIST.  However,  only  the  first  two  quantities  are  used  for  the  actual 
tracking  calculations.  The  track  is  accurate  to  within  a few  arc  seconds  and 
the  star  is  always  observed  in  the  center  of  the  field  of  view  of  the  main 
telescope,  as  long  as  current  RA  and  DEC  values  from  the  catalog  are  supplied. 

The  SATTRACK  routine  is  a passive  satellite  tracking  program  which 
requires  the  two-card  parameter  set  provided  by  NORAD  as  an  input.  It  uses 
the  Simplified  General  Perturbation  (SGP)  theory  to  calculate  the  instantaneous 
satellite  position.  The  satellite  offsets  from  the  telescope  center  vary 
from  satellite  to  satellite  and  depend  on  the  accuracy  of  the  card  parameters. 
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The  mount  control  commands  described  in  Section  2.4.  of  the  User's  Manual  can 
then  be  used  to  introduce  offsets  in  the  calculated  mount  position  to  bring 
the  satellite  to  the  center  of  the  field  of  view.  Sometimes  the  offset  is 
large  enough  for  the  satellite  not  to  be  visible  in  the  main  scope,  which  has 
a field  of  view  of  only  5 arc  minutes,  but  it  is  easily  visible  in  the  viewing 
telescope  which  has  a 1°  field  of  view.  None  of  the  satellites  that  were 
tracked  caused  the  mount  to  slew;  i.e.,  the  speed  of  the  satellite  was  well 
within  the  capability  of  the  mount. 

The  KALTRACK  routine  is  the  active  tracking  program  which  requires  a 
number  of  approximate  azimuth  and  elevation  angles  of  the  satellite  to  be 
supplied  manually  at  sequential  time  instants.  The  program  uses  these  input 
measurements  to  update  the  satellite  parameters.  In  a follow-up  contractual 
work,  the  necessity  to  input  manually  will  be  replaced  by  hardware.  Since 
updating  of  parameters  is  possible  only  after  taking  repeated  measurements, 
the  initial  parameter  set  with  which  the  program  is  started  must  be  accurate 
enough  for  the  satellite  to  be  visible  in  the  telescope. 

6.2.  SYSTEM  DESCRIPTION 

This  section  describes  the  system  with  the  aid  of  flowcharts  at  different 
levels  of  complexity. 

6.2.1.  Overall  System 

Figure  6-1  is  a flowchart  of  the  operations  required  to  bring  up  the 
Real-time  Disk  Operating  System  (RDOS)  on  the  NOVA  800  computer.  A further 
detailed  description  is  available  in  the  User's  Manual  to  complement  the 
flowchart. 


Depending  on  the  tracking  operation  to  be  performed,  the  user  can  type 
STARTRACK,  SATTRACK  or  KALTRACK  to  read  the  appropriate  program  into  the 
core,  Figure  6-2.  The  program  in  core,  after  getting  the  appropriate  para- 
meters and  commands,  generates  azimuth  and  elevation  angles  as  functions  of 
time  and  transmits  to  the  mount  through  the  mount-computer  hardware  interface. 


I 

f 


6.2.2.  File  Maintenance 

Each  of  the  three  tracking  systems  (STARTRACK,  SATTRACK  and  KALTRACK) 
maintains  a data  file  on  the  objects  tracked  frequently.  The  command  LIST 
(2. 3. 2. 2.)  may  be  used  to  retrieve  the  STARLIST  or  SATLIST  as  the  case  may 
be.  Modification  of  the  list  resident  in  the  core  is  possible  by  reading  in 
cards.  Using  the  RDOS  command  CARDINPUT),  for  SATTRACK  or  KALCARD  for 
KALTRACK,  a reorganized  and  updated  listing  of  the  list  is  printed  on  the 
line  printer  for  either  command  (Figure  6-3,  6-4).  To  obtain  a printout  of 
the  KALLIST,  it  is  necessary  to  use  the  command  KALCARD V with  maximum  number 
of  cards  as  zero.  KALLIST  cannot  be  obtained  while  the  tracking  is  in  pro- 
gress, as  opposed  to  STARLIST  and  SATLIST  which  can  be  obtained  during  track 
by  using  the  command  LIST  . The  STARLIST  can  accomodate  100  star  para- 
meters. Each  set  of  parameters  for  a star  consists  of  the  catalog  number  of 
the  star,  brightness  magnitude,  right  ascension  in  hours,  minutes,  seconds 
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End  of 
File  Return 


Read  a set  of  card  parameters  from  the  list 
on  the  disk.  If  the  end  of  file  mark  is 
encountered,  take  end  of  file  return. 


Figure  6-3  Flow  Chart  for  Card  Input 
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Figure  6-3  (Continued) 
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Figure  6-4  Flowchart  for  Kalcard 
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and  declination  in  degrees,  minutes,  seconds.  The  SATLIST  can  accomodate 
parameter  sets  for  50  satellites,  12  parameters  per  satellite.  The  set  for 
one  satellite  consists  of  24  parameters  in  Kaltrack.  Kaltrack  can  read  in 
any  number  of  satellites.  The  file  maintenance  procedure  for  the  three 
programs  is  represented  in  Figures  6-5,  6-6  and  6-7. 

6.2.3.  Tracking  Systems 

During  any  of  the  three  tracking  operations,  two  programs  are  being 
executed.  The  main  program,  which  calculates  the  tracking  parameters,  has 
the  higher  priority  and  is  being  executed  as  often  as  possible.  The  second 
program,  with  a lower  priority  rating,  is  a Command  Line  Interpreter  which 
gets^gxecuted  only  in  available  spare  time.  The  clock  supplies  a pulse  every 
1/60  of  a second.  Each  star-tracking  calculation^akes  less  than  one  clock 
pulse  period  and  the  time  left  over  out  of  the  1/60T  second  is  used  to 
execute  the  CLI.  A satellite-tracking  calculation  takes  slightly  over  five 
clock  pulse  periods.  The  remaining  time  in  the  sixth  pulse  period  is  used 
for  executing  the  CLI.  The  CLI  can  accept  four-letter  commands  from  the 
teletype  and  set  or  clear  flags  accordingly.  (Figure  6-8).  These  flags 
determine  the  course  of  the  main  program.  All  acceptable  commands  have  been 
described  in  the  User's  Manual. 

6. 2. 3.1.  STARTRACK 

The  program  for  tracking  stars  has  been  described  in  PAR  Report  #75-32 
on  the  first  phase  of  this  project  (Contract  #F30602-75-C-0144) . A short 
description  is  repeated  here  for  the  sake  of  completeness.  Figure  6-9  is  a 
functional  flowchart  of  the  star-tracking  program.  The  program  reads  in  the 
right  ascension  (RA)  and  declination  (DEC)  angles  of  the  star  to  be  tracked 
and  defines  a unit  vector  pointing  in  that  direction  in  the  earth-centered 
inertial  (ECI)  coordinates.  Coordinate  transformation  from  ECI  axes  to  site 
topocentric  axes  is  performed  and  corresponding  azimuth  and  elevation  angles 
are  calculated.  This  transformation  makes  use  of  the  clock  to  compute  the 
sidereal  time  and  get  the  current  site  vector.  Azimuth  and  elevation  offsets 
are  then  transferred  to  the  mount  servos.  Then,  RDOS  tasks  the  CLI  to  check 
if  any  commands  came  in  via  the  teletype.  These  operations  are  repeated  at 
almost  60  Hz. 

The  STARTRACK  is  an  accurate  program  and  it  keeps  the  telescope  locked 
on  to  the  star. 

6. 2. 3. 2.  SATTRACK 

SATTRACK  is  essentially  a satellite  position  prediction  routine  based  on 
a Simplified  Generalized  Perturbation  (SGP)  theory  described  in  Section  3. 

The  SATLIST  provides  the  orbital  parameters  at  some  epoch.  The  prediction 
routine  calculates  the  updated  position  in  the  orbit  by  applying  the  secular, 
long-period  and  short -period  perturbations  and  then  applies  the  transforma- 
tions to  — first  the  ECI  coordinates,  and  then  to  the  local  topocentric 
coordinates  (Figure  6-10). 


Figure  6-5  File  Maintenance  for  STARTRACK  Routine 
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Figure  6-6  File  Maintenance  for  SATTRACK  Routine  "SATLIST" 
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Figure  6-7  File  Maintenance  for  Kalman  Filter  Routine  "KALLIST" 


Using  kalcard  the  list  on 
the  disk  can  be  updated  by 
deleting  specific  satellite! 
from  the  disk  by  list  #. 
Also  cards  can  be  read  in. 
For  a duplicate  satellite 
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list  as  necessary. | 
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done  filtering  he 
can  write  the  new 
i parameters  out  to 
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given  list  position. 
If  the  position  he 
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disk,  the  param- 
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after  the  last  list 
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cannot  delete  a 
satellite  from  the 
list  in  the  track- 
ing routine. 
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Figure  6-8 


Operation  of  Tracking  System?;  under  RDOS 
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\ adjusted. 
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Positions  are 
output  to  mount 
driver. 
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This  side  is  of  lowest 
priority  and  is  executed 
only  in  the  spare  time  of 
the  track  routine. 


(A  record  of  the  track  is 
kept  and  command  line 
interpreter  is  tasked. 


The  other  side  is  of 
highest  priority  and  is 
continuously  being 
executed. 


Track  is  cued 
to  stop  and  start 
at  correct  time. 

One  clock  of  delay 
is  added^to  system. 
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Figure  6-9  Flowchart  for  Task  1 
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6. 2. 3. 3.  KALTRACK 

When  it  is  possible  to  measure  the  actual  azimuth  and  elevation  of  the 
satellite  as  a function  of  time,  KALTRACK  can  be  used  to  track  the  satellite 
and  update  the  orbital  parameter,  simultaneously.  A NOVA-TV  interface  has 
been  designed  and  fabricated  for  the  specific  purpose  of  measuring  the 
actual  satellite  differential  position  after  it  is  acquired  in  the  field  of 
view.  This  hardware  is  described  in  Section  7 of  this  report. 

According  to  the  theory  described  in  Section  5,  the  satellite  parameters 
may  be  updated  using  the  measured  azimuth-elevation  angles  of  the  satellite 
using  the  Kalman  filter  technique  (Figure  6-11). 

6.3.  HARDWARE 

6.3.1.  Computer-Mount  Interface 


Only  the  computer-mount  interface  designed  by  PAR  will  be  described  in 
short  for  the  sake  of  completeness.  A detailed  design  description  was  sub- 
mitted as  an  Interim  Report.  Figure  6-12  schematically  shows  the  hardware- 
software  interface  chart.  Supplying  the  new  look  angles  to  the  mount,  and 
reading  the  current  mount  position,  are  the  two  functions  performed  through 
the  interface.  Error  codes  are  also  generated  which  indicate  whether  or  not 
the  mount  is  slewing  or  at  an  endstop. 

Figure  6-13  is  a schematic  of  the  data  link  between  the  CPU  and  the 
mount  for  relaying  the  new  look  angles  to  the  mount.  Figure  6-14  shows  the 
mount  to  CPU  interface  for  reading  the  current  mount  position  and  6-15  the 
clock  CPU  interface. 


6.3.2.  Computer-TV  Interface 

A computer-TV  interface  was  fabricated  to  measure  the  actual  satellite 
position.  When  the  satellite  is  acquired  on  the  TV  screen,  cross-hairs  on  the 
screen  can  be  moved  to  its  observed  position,  and  a push  button  switch 
supplies  these  azimuth-elevation  offsets  to  the  Kalman  filter  as  inputs.  The 
filter  applies  a calculated  correction.  This  procedure  can  be  repeated  by 
moving  the  cross-hairs  to  the  new  satellite  position  The  hardware  and  the 
operation  has  been  described  in  Section  7. 


6.4. 


SUBROUTINE  AND  PROGRAM  FUNCTION  AND  INDEX 


6.4.1. 


STARTRACK 


The  routines  described  in  this  subsection  are  exclusive  to  the  STARTRACK 
program.  STARTRACK  is  the  main  driver  which  calls  the  appropriate  routine  to 
implement  the  track. 
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Filter? 


Compute  Sensitivities  of' 
Orbital  Parameters  with 
Respect  to  AZ  EL  Error 


Get  the  AZ  EL 
Error 


Figure  6-11  Kalman  Filter  Routine 


(Multiply  the  AZ  EL  Error  ^ 
Times  the  Sensitivities  to 
get  Corrections  to  the 
Vorbital  Parameters  ) 


Add  the  Correc-  \ 
tions  to  the  j 
Orbital  Parameters 


yftre  ' 
'we  to 


'Update  PK? 


Figure  6-11  (Continued) 
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,/CPU  Loads  Four' 
Words  to  Inter- 
i face  Buffers 

V V 


CPU  Gives  ' 

Transfer  Commanc 


EL  Low  Word 


The  multiplexer  has  the  AZ  and  EL  positions  as  its 
inputs  at  all  times.  Upon  CPU  command,  the  data  is 
strobed  onto  the  mount  CPU  cable.  When  AZ  data  is  on 
the  cable,  the  AZ  load  is  high.  The  same  for  EL. 


CPU  to 

Mount  Date 

> Cable 


18  Bit  Mount  Position 


AZ  Mount 
Position 


EL  Mount 
Position 


Mount  End  of  Data  Link 


Figure  6-13  CPU  End  of  Data  Link 


EL  poad 


Figure  6-14  Mount  to  CPU  Interface  to  Input  Mount  Position 


AZ  Parallel 
Position  Buffer 


EL  Parallel 

Position 

Buffer 


18  Bits 


AZ  Shift 
Register 


EL  Shift 
Register 


Mount  Side  of  Interface 
CPU  Side  of  Interface 


jAZ  Shift 
■Register 


EL  Shift 
Register 


'CPU  Request 
Path  Transfer 


Controller  First  Sends  a 
Pulse  to  Load  the  Mount 
Shift  Register 


Controller  Sends  18  Shift 
Pulses  to  Serially  Load  the 
Position  Buffers  Out  the 
Interface  Board 


ft PU  Inputs  the 
Data  as  Four 
, Words 
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Figure  6-15  Clock  CPU  Interface 
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6. 4. 1.1.  C0MND1 


This  is  the  command  line  interpreter  for  the  STARTRACK  program  which  has 
no  real  input  or  output,  but  serves  to  modify  parameters  and  flags  of  other 
routines,  according  to  the  operator's  commands. 

6. 4. 1.2.  PAR3 

This  routine  initializes  the  stellar  tracking  routine  by  supplying 
output  parameters  in  the  common  areas. 

COMMON/ SDRIVE/  # 

[RX,  RY,  RZ] 

/II/ 

IYR 
I DAY 
/ 12/ 

LAT 

COSLAT,  SINLAT 
ST  IKE 
OSTIME 
OTSINCE 
SINS,  COSS 
VXY,  VZ 

6.4.2.  SATTRACK4 

This  routine  links  all  the  appropriate  routines  to  do  a passive  satellite 
track.  All  the  subroutines  exclusive  to  this  system  have  been  described  in 
this  subsection. 

6. 4. 2.1.  COMND 

This  is  the  SATTRACK  analog  of  C0MND1  and  is  the  command  line  interpreter 
for  the  passive  satellite  tracker  routine.  It  has  no  real  inputs  or  outputs 
except  that  it  modifies  parameters  and  flags  in  response  to  operator  commands 
to  control  the  course  of  other  subroutines. 

6. 4. 2. 2.  ORBIT 

This  routine  performs  the  transformation  of  parameters  from  the  satellite 
orbital  coordinate  system  to  the  earth-centered  inertial  system  (ECI)  parameters. 


Initial  star  position 
Year 

Number  of  day  in  the  year 
Site  latitude 

Latitude  cosine  and  sine  functions 
Current  sidereal  time 
Epoch  time 
Time  since  epoch 

Sine  and  Cosine  of  the  sidereal  time 
Related  to  flatness  of  earth 


Input  Parameters 

TRUEU  A dummy  variable  (true  anomaly) 

/05/ 

rmAG  Magnitude  of  satellite  vector  in  ECI  axes 

^ The  name  within  the  slashes  /SDRIVE/  is  the  name  of  the  common  block. 
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SNODE,  CNODE 
SINI,  COSI 
RMGDT 
RVDT 
1171 
COSLAT 
SINS,  COSS 
CVXY , CVZ 
Output  Parameters 
/SDRIVE/ 

RX,  RY,  RZ 

6. 4. 2. 3.  PAR 5 

This  routine  performs  a similar  function  in  SATTRACK  as  does  PAR3  in 
STARTRACK.  It  has  some  I/O  with  the  operator  which  involves  initializing  the 
program  as  described  in  the  User’s  Manual  under  2. 2. 2. 4. 

6. 4. 2. 4.  SEMI  (Eccentricity,  Mean  Motion,  Inclination) 

This  function  computes  the  semi-major  axis  of  the  orbit  given  the  eccen- 
tricity, mean  motion  and  the  inclination.  It  is  also  used  in  KALTRACK. 


Sine  and  cosine  of  (RA  of  the  ascending  node) 
Sine  and  cosine  of  i (inclination) 

Magnitude  of  rate  of  change  of  RMAG 
Radial  component  of  velocity 

COS  (Latitude) 

Sine  and  cosine  of  the  sidereal  time 
Constants  due  to  oblateness  of  earth 


X,Y,  Z in  the  ECI  system 


Input  Parameters  as  arguments  are  the  dummy  variables  representing 
eccentricity,  mean  motion  and  inclination  as  shown  above. 


Output  Function 
SEMI 

6. 4. 2. 5.  TRANS  (AZ,  EL) 


Semi-major  axis 


This  routine  serves  to  perform  the  transformation  from  the  ECI  to  site 
topocentric  axes  and  then  to  calculate  the  azimuth  and  elevation  angles. 


Input  Paramters 

/SDRIVE/ 

RX,  RY,  RZ 
/ 1 2 / 

COSLAT,  SINLAT 
SINS,  COSS 
CVXY,  CVZ 
LAT 

Output  Parameters 
AZ,  EL 


Position  of  object  in  ECI  coordinate  system 

Sine  and  cosine  of  latitude 
Sine  and  cosine  of  sidereal  time 
Constant  for  flatness  of  earth 
Latitude 


AZ,  EL  positions  of  objects 
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6. 4. 2. 6.  VARI  (TRUE,  TSINCE) 


Using  the  classical  orbital  parameters  and  the  time  since  epoch  (TSINCE), 
this  routine  calculates  the  satellite  position  in  the  orbital  coordinate 
system.  Common  area  /II/  is  the  input  to  this  routine  and  the  output  para- 
meters are  loaded  into  the  common  area  /01/,  /02/,  /05/,  /TEST/ 

/II/ 

EO,  EPOCH,  10,  MO,  NDOTO , ND0T6 , NO,  NODEO,  OMEGO,  REVO 
/SDRIVE/ 

RX,  RY,  RZ 

/01/  ' 

IFLG 

AM,  EM,  IM,  NODEM,  0MEGM,J>  Mean  Values  of  parameters 
NDOTM,  REVM,  LM,  NM 
/02/ 

ELONG,  LLONG, 

EXLNG , OMEGL 
TEMP,  TEMPL 
/TEST/ 

IS,  NODES 
/05/ 

RMAG,  SNODE,  CNODE, 

RMGDT,  RVDT 

6. 4. 2. 7.  CARD 

This  routine  reads  the  NORAD  two-card  data  set  into  the  core  and  into 
the  following  common  areas  and  sets  error  flag  (FLAGE)  if  there  is  an  error. 

/II/  EPOCH,  YR,  REVO,  NO,  10,  EO,  MO,  NODEO,  OMEGO, 

NDOTO,  ND0T6 
/CARDC/  SATNO,  FLAGE 

6. 4. 2. 8.  CARDINPUT 

This  is  a main  program  which  updates  the  satellite  list  by  introducing 
the  satellite  data  read  from  the  cards  in  their  proper  order  in  the  SATLIST. 

If  the  data  of  any  satellite  being  read  in  already  exists  in  SATLIST,  then 
the  old  data  is  overwritten  with  the  new  parameter  set  (2. 2. 2. 2.).  The 
variables  used  are  all  those  listed  under  CARD. 

6.4.3.  KALTRACK 

This  program  uses  Linearized  Kalman  Filter  to  perform  an  active  satellite 
track,  with  the  operator  supplying  azimuth  and  elevation  angle  measurements 
of  the  observed  satellite.  The  observations  are  used  to  form  an  error  signal 
which  is  used  to  apply  a correction  to  the  orbital  elements.  When  there  are 
no  observations  forthcoming,  KALTRACK  uses  a variation  of  SATTRACK  to  continue 
tracking  the  satellite  by  prediction  alone.  KALTRACK  is  the  driver  program 
to  perform  the  tracking  function  and  has  no  inputs  and  outputs  of  its  own. 

Most  variables  occur  in  this  program. 


^ Long-period  perturbed  values 

J 


Short-period  perturbed  values 


SINI , COSI 
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6. 4. 3.1.  EXANM  (Mean  Anomaly,  Eccentricity) 

This  function  is  used  in  both  KALTRACK  and  SATTRACK4  to  calculate  the 
eccentric  anomaly  using  the  eccentricity  and  mean  anomaly.  Kepler's  equation 
is  solved  by  employing  a simple  iteration  technique. 

Input  Parameters  are  the  arguments  Mean  Anomaly  and  Eccentricity  as 
shown  above. 

Output  Function 

EXANM  Eccentric  anomaly 

6.4. 3. 2.  KAL  (TSINCE,  A0,  E0) 

This  is  the  main  Kalman  filter  routine  which  updates  the  orbital  para- 
meters based  on  the  azimuth  and  elevation  measurements  and  transmits  the 
predicted  az-el  angles  to  the  mount  for  tracking  purposes.  It  detects  any 
mount  or  time  offsets  from  the  system. 

Input  Parameters 

TSINCE 
/II/ 

/ 1 2/ 

/KALF/ 

Output  Parameters 

/II/ 

/OFSET/ 

A0 , E0 

/01/ 

6.4. 3. 3.  KALCOMND 


EPOCH,  YR,  REVO,  NO,  10,  EO,  MO,  NODEO^,  Card 
OMEGO,  NDOTO,  ND0T6,  IYR,  DAY  parameters 

COSLAT,  SINLAT,  STIME,  OTSINCE,  VXY,  VZ'y  Earth's 
SINS,  COSS,  LAT  / position 

XHATK,  PK,  QX,  CRK1,  NCDOT,  NCK,  DELTK,  KALFLAG 


EPOCH,  YR,  NO,  10,  EO,  MO  NODEO,  OMEGO 
OAZ,  OEL,  OT  System  offset 
AZ,  EL  Look  angles 
EM 


Card 

parameters 


This  routine  is  the  command  line  interpreter  for  KALTRACK.  It  has  no 
real  inputs  or  outputs,  but  it  modifies  the  parameters  and  flags  of  other 
routines . 


6.4. 3.4.  KALPAR 


This  routine  intializes  the  Kalman  filter  tracking  system  with  the 
following  output  parameter. 


/II/ 
/ 1 2/ 


EPOCH,  YR,  REVO,  NO,  10,  EO,  MO,  NODEO, 
OMEGO,  NDOTO,  ND0T6 , IYR,  DAY 
COSLAT,  SINLAT,  STIME,  LAT,  OSTIME, 
OTSINCE,  VXY,  VZ,  SINS,  COSS 
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6.4. 3. 5.  KALP 


(TRUIU,  TSINCE,  XBARK,  R0WTK1,  AZ0 , EL0) 


This  routine  performs  the  prediction  part  of  the  filter.  Using  the  most 
recent  orbital  parameter  estimate  it  predicts  the  state  transition  matrix  and 
the  az-el  angles. 


Input  Parameters 


TSINCE 

/II/ 


/ 12  / 


TIME  science  epoch  . 

EPOCH,  YR,  REVO,  NO,  10,  E0,  M0^>  Card 
N0DE0,  0MEG0,  ND0T0,  ND0T6  parameters 

COSLAT,  SINLAT,  LAT,  CVXY,  CVZ,  SINS,^>  Earth's 
COSS  position 


Output  Variables 


/01/ 

/02/ 

/05/ 

/TEST/ 

/SDRIVE/ 


TRUEU,  XBARK(6),  R0WTK1,  AZ0,  EL0 
6. 4. 3. 6.  KALCARD 


OMEGM,  NODEM,  NM,  EM,  REVM,  IFLG,  IM,  LM,  NDOTM 
OMEGL,  TEMP,  TEMPL,  ELONG,  EXLNG,  LLONG 
RVDT,  CNODE , COSI,  SINI,  RMAG,  RMGDT 
IS,  NODES 
RX,  RY,  RZ 


This  is  a main  program  which  updates  the  KALLIST  according  to  the  new 
card  parameter  set  fed  in. 


6. 4. 3. 7.  OBSER  (AZ,  EL) 

This  routine  passes  the  observed  angles  to  KAL. 

6. 4. 3. 8.  PKQX  (CARD1,  CARD2,  PK,  QX) 


This  routine  supplies  the  initial  values  for  the  initial  error  covariance 
PK  and  the  system  noise  covariance  QX  based  on  the  new  and  old  card  parameters. 


Input  Parameters 


CARD1 
CARD  2 


Old  set  of  card  parameters 
New  set  of  card  parameters 


Output  Parameters 


PK 

QX 


Error  covariance  of  orbital  parameters 
System  noise  covariance  matrix 


6.4.4.  General  Purpose  and  Common  Routines 


This  subsection  describes  the  routines  which  are  not  exclusive  to  any 
one  of  the  above  three  systems  and  perform  frequently  required  general- 
purpose  functions. 
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6. 4. 4.1.  Matrix  Operations  I 

6. 4. 4. 1.1.  ABAT  (A,  B,  C,  I,  J,  D)  1 

This  provides  the  following  matrix  operation  1 

C = ABAT 

where  A and  B are  (IxJ)  and  (JxJ)  input  matrices  and  C is  the  (Ixl)  output 
matrix.  D is  an  (IxJ)  matrix  used  for  intermediate  result. 

6. 4. 4. 1.2.  ABT  (A,  B,  C,  I,  J,  K) 

This  performs  the  matrix  operation 

T 

C = AB 

where  A and  B are  (IxJ)  and  (JxK)  input  matrices  and  C is  the  (IxK)  output 
matrix. 

6. 4. 4. 1.3.  MATMUL  (A,  B,  C,  I,  J,  K) 

As  the  name  implies,  this  routine  performs  the  matrix  product,  C=AB 
where  C,  A and  B are  (I,K),  (I,J)  and  (J,K)  matrices  respectively. 

6. 4. 4. 1.4.  MATSUB  (A,  B,  C,  I,  J) 

This  routine  outputs  the  difference  of  two  (I,J)  matrices  A and  B as  C. 

6. 4. 4. 1.5.  MATADD  (A,  B,  C,  I,  J) 

This  routine  outputs  the  sum  of  two  matrices  A and  B as  C,  all  being  of 
dimensions  (I,J). 

6. 4. 4. 2.  TRIGONOMETRIC  OPERATIONS 

6. 4. 4. 2.1.  M0D2P  (X) 

This  routine  takes  the  M0D(2  tt  ) of  a number  X. 

6. 4. 4. 2. 2.  AKTAN  (XX,  YY) 

This  function  takes  the  arc  tangent  when  the  x and  y components  are 
supplied  to  it  as  dummy  variables.  The  output  angle  is  between  0 and  2 it. 

6. 4. 4. 3.  CLOCK  ROUTINES 

All  routines  having  anything  to  do  with  time  have  been  categorized  under 
this  heading. 
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6.4.4. 3.1.  TIME  (TSINCE) 


This  routine  calculates  the  time  since  the  last  update  of  the  orbital 
parameters  and  the  sidereal  time  using  the  current  time  of  day  and  values  of 
those  parameters  at  the  beginning  of  the  day. 

Input  Parameters 

/OFSET/  OT  (System  time  offset) 

/I2/  OSTIME,  OTSINCE  (Sidereal  time  and  tsince  at 

24:00  the  day  before.) 

Output  Parameters 


TSINCE 

/I2/ 


Time  since  epoch 

STIME,  SINS,  COSS  (Sidereal  time  and  its 
sine  and  cosine.) 


6. 4. 4. 3. 2.  CLCK  (MS) 


CLCK  returns  the  time  of  day  in  milliseconds. 

6. 4. 4. 3. 3.  CLCK1  (W0RD1,  W0RD2) 

CLCK1  returns  the  time  of  day  as  two  words:  i.e.,  tens  of  seconds  and 

milliseconds.  Output:  first  two  locations  on  the  stack. 

6. 4. 4. 4.  SYSTEM  ROUTINES 

6. 4. 4. 4.1.  FDELY  (X) 

This  is  an  RDOS  routine  which  can  be  employed  to  delay  the  advance  of 
the  program  by  a desired  number  of  clock  pulses. 

6. 4. 4. 4. 2.  FTASK  (SUBROUTINEJJAME,  $ERROR_RETURN,  PRIORITY) 

This  is  another  RDOS  routine  which  sets  up  a task  control  block  for  a 
new  task. 

6.4.4. 5.  MOUNT-COMPUTER  INTERFACE  ROUTINES 

6. 4. 4. 5.1.  MNTS 


This  routine  starts  the  mount-computer  communications  by  synchronizing 
the  two  devices. 

6.4.4. 5.2.  MOUNT  (ERR) 


J 

i 

« 

j 

1 

J 

h 

H 
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t 
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This  routine  converts  the  input  angle  from  radians  to  revolutions  and 
outputs  the  result  to  the  assembly  language  driver  routine. 
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Input  Parameters 
/MNT/ 


AZ0,  EL0,  (AZ,  EL  look  angles) 
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Output  Parameters 

/MNT/  TIME2 , AZ1,  ELI  (AZ,  EL  position  and  time 

in  days  of  position  read) 

6.4.4. 5. 3-;  MNT1  (AZ,  EL) 

This  is  an  assembly  language  driver  to  output  the  angles  to  the  mount. 
Input  parameters  are  the  first  two  locations  on  the  stack. 

6. 4. 4. 5. 4.  MNT 2 (AZ,  EL,  AZER,  ELER,  FLAG2 , FLAG3) 

This  also  is  an  assembly  language  driver  which  accepts  input  from  the 
mount  in  the  first  six  locations  on  the  stack. 

Output  Parameter 

AZ,  EL  AZ,  EL  position  in  revolutions 

AZER,  ELER  AZ,  EL  hardware  error  positions 

FLAG2,  FLAG3  Video  interlace  control  bits 

6.4.4. 5. 5.  ERROR  (ERR) 

(ERR)  Input  error  code 

This  routine  decodes  any  errors  incurred  during  the  positioning  of  the 
mount . 

6. 4. 4. 6.  SEARCH  (Input  Azimuth,  Input  Elevation,  Output  Azimuth,  Output 
Elevation,  Mode  Flag) 

This  routine  implements  a spiral  search  in  the  line  of  sight  coordinate 
system. 

MODE  = 0 if  search  calculation  is  to  be  done 

= 1 is  search  is  to  be  setup  6 controlled 
= 2 if  search  is  to  be  controlled 

Input  Parameters 

/SRCH/  6 /OFSET/  OAZ,  OEL,  TIME 

Output  Parameters 

/OFSET/  OAZ,  OEL  system  AZ,  EL  offsets 
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6.4.5. 

Alphabetical  List  of 

Subroutines 

Number 

Name 

Section  Number 

Page  Number 

1 

ABAT 

6. 4. 4. 1.1. 

2 

ABT 

6. 4. 4. 1.2. 

6-31 

3 

AKTAN 

6. 4. 4. 2. 2. 

6-31 

4 

CARD 

6. 4. 2. 7. 

6-28 

5 

CARD INPUT 

6. 4. 2. 8. 

6-28 

6 

CLCK 

6.4.4. 3.2. 

6-32 

7 

CLCK1 

6. 4. 4. 3. 3. 

6-32 

8 

COMND 

6. 4. 2.1. 

6-26 

9 

C0MND1 

6. 4. 1.1. 

6-26 

10 

ERROR 

6. 4.4. 5. 5. 

6-33 

11 

EXANM 

6. 4. 3.1. 

6-28 

12 

FDELY 

6. 4. 4. 4.1. 

6-32 

13 

FTASK 

6. 4. 4. 4. 2. 

6-32 

14 

KAL 

6. 4. 3. 2. 

6-29 

15 

KALCARD 

6. 4. 3. 6. 

6-30 

16 

KALCOMND 

6. 4. 3. 3. 

6-29 

17 

KALP 

6. 4 . 3.  5 . 

6-30 

18 

KALPAR 

6. 4. 3. 4. 

6-29 

19 

KALTRACK 

6.4.3. 

6-28 

20 

MATADD 

6. 4. 4. 1.5. 

6-31 

21 

MATMUL 

6. 4. 4. 1.3. 

6-31 

22 

MATSUB 

6. 4. 4. 1.4. 

6-31 

23 

MOD2P 

6. 4. 4. 2.1. 

6-31 

24 

MOUNT 

6.4.4. 5.2. 

6-32 

25 

MNTS 

6. 4. 4. 5.1. 

6-32 

26 

MNT1 

6. 4. 4. 5. 3. 

6-33 

27 

MNT2 

6.4.4. 5.4. 

6-33 

28 

OBSER 

6. 4. 3. 7. 

6-30 

29 

ORBIT 

6. 4. 2. 2. 

6-26 

30 

PAR3 

6. 4. 1.2. 

6-26 

31 

PAR5 

6.4. 2. 3. 

6-27 

32 

PKQX 

6.4. 3. 8. 

6-30 

33 

SATTRACK4 

6.4.2. 

6-26 

34 

SEARCH 

6. 4. 4. 6. 

6-33 

35 

SEMI 

6.4. 2.4. 

6-27 

36 

STARTRACK 

6.4.1. 

6-19 

37 

TIME 

6.4.4. 3.1. 

6-32 

38 

TRANS 

6. 4. 2. 5. 

6-27 

39 

VARI 

6. 4. 2. 6. 

6-27 

6.4.6. 

Alphabetical  List  of 

I/O  Variables 

Math  Equivalent 

Occurs  in  Program 

A 

a - semi  major  axis 

KALP,  VARI 

AM 

a -mean 

KALP,  VARI 

AO 

a-original 

KALP,  VARI,  KAL 

AP 

AZ  predicted  angle 

SEARCH 
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A0 

Azimuth  pointing  angle 
without  any  offsets 

KAL 

AZ 

Azimuth 

Most  programs 

AM 

Mean  value  of  semi-major 
axis 

VARI,  KALP 

AZ0 

Azimuth  output  angle 

KALTRACK,  STARTRACK, 
MOUNT,  SATTRACK,  KALP 

AZ1 

Azimuth  input  angle 

KALTRACK,  STARTRACK, 
MOUNT,  SATTRACK 

AZER 

Azimuth  error  from  hardware 
interface 

MNT2 

CARD1 

First  set  of  card  parameters 

PKQX 

CARD  2 

Second  set  of  card  parameters 

PKQX 

CNODE 

Cos  (node) 

KALP,  ORBIT,  VARI 

COSI 

Cos  (inclination) 

KALP,  ORBIT,  VARI 

COSLAT 

Cos  (latitude) 

KAL,  TRANS,  KALP,  ORBIT 
PAR3 , PAR5 , KALPAR 

COSS 

Cos  (Sidereal  time) 

TIME,  TRANS,  KALP,  ORBIT 
PAR3 , PAR5 , KALPAR 

CRK1 

Covariance  of  AZ,  EL  inputs 

KAL 

CVXY 

Constant  due  to  oblateness 
of  earth 

KALP,  TRANS,  ORBIT 

CVZ 

Constant  due  to  oblateness 
of  earth 

KALP,  TRANS,  ORBIT 

DAY  (or  IDAY) 

Day  of  year 

KALPAR,  PAR3 , PARS,  KALPAR 
Tracking  routines,  KAL 

DELTK 

Time  since  last  Kalman  filter 

KAL 

E0 

Elevation  pointing  angle 
without  any  offsets 

KAL 

EL 

Elevation 

Most  programs 

EL0 

Elevation  output  angle 

Main  tracking  routine 
and  MOUNT,  KALP 

ELER 

Elevation  error  from 
hardware  interface 

MNT2 

ELONG 

Long  periodic  term  of 
eccentricity 

KALP,  VARI 

ELI 

Elevation  input  angle 

Main  Tracking  routines 
and  MOUNT 

EM 

Mean  eccentric  anomaly 

KALP,  KAL,  VARI 

EO 

Original  Eccentricity 

KALPAR,  KALP,  KAL, 
VARI,  CARD,  CARD INPUT 

EPOCH 

Time  of  original  values 

KALPAR,  CARD,  KALP, 
KAL,  VARI,  CARDINPUT 

ERR 

Error  code  for  mount 
positioning 

MOUNT,  ERROR 

EXLNG 

Long  periodic  term  for 
eccentric  anomaly 

KALP,  VARI 

FLAG2,  FLAG3 

Flags  from  hardware  interface 

MNT2 

FLAGE 

Error  code  for  subroutine 
"CARD" 

CARD,  KALCARD,  CARDINPUT 

IDAY 

Sa.ne  as  DAY 

IFLG 

Used  to  short  cycle  VARI 

VARI,  KALP 
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i 

x 

t 


i 


IM 

Mean  inclination 

10 

Original  inclination 

IS 

Short  periodic  term  for 
inclination 

IYR  (or  YR) 

Year 

KALFLAG 

Flag  to  start  Kalman  filter 

LAT 

Latitude 

LLONG 

Orbital  longitude 

LM 

Mean  orbital  longitude 

MO 

Mean  anomaly 

MODE 

Mode  flag  for  search  routine 

NCDOT 

Corrected  (estimated)  NDOTO 

NCK 

Unused  variable 

NDOTM 

Mean  value  of  NDOTO 

NDOTO 

First  derivative  of  mean 
motion 

ND0T6 

Second  derivative  of  mean 
motion 

NM 

Mean  value  of  mean  motion 

NO 

Mean  motion  at  epoch 

NODEM 

Mean  value  of  NODE 

NODEO 

Right  ascension  of  ascending 
node 

NODES 

Short  periodic  term  of  NODE 

OAZ 

Azimuth  offset  for  system 

OEL 

Elevation  offset  for  system 

OMEGL 

Long  periodic  term  for  Omega 

OMEGM 

Mean  periodic  term  for  Omega 

OMEGO 

Original  Omega 

OSTIME 

Sidereal  time  at  the 
beginning  of  the  day 

OT 

Time  offset  for  system 

OTSINCE 

Time  since  epoch  at  beginning 
of  the  day 

PK 

Confidence  matrix  for  the 
system 

QX 

Error  covariance  matrix  of 
system  noise 

REVM 

Mean  number  of  revolutions 
from  launch 

REVO 

Number  of  revolutions  since 
launch  to  epoch 
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VARI,  KALP 
KAL,  VAP.I , KALP, 

CARD,  CARDINPUT,  KALPAR 
VARI,  KALP 

CARDINPUT,  CARD,  main 
tracking  routines,  PAR3, 
PAR 5,  KALPAR 
KAL 

KALP,  KALPAR,  TRANS, 

KAL,  PAR3 , PAR5,  KALPAR 
VARI,  KALP 
VARI,  KALP 
KALPAR,  VARI,  KALP, 

CARD,  CARDINPUT,  KAL 

SEARCH 

KAL 

KAL 

VARI,  KALP 
VARI,  KALP,  CARD, 
CARDINPUT,  KAL,  KALPAR 
VARI,  KALP,  CARD, 
CARDINPUT,  KAL,  KALPAR 
VARI,  KALP 

KAL,  VARI,  KALP,  CARD, 
CARDINPUT,  KALPAR 
VARI,  KALP 
VARI,  KALP,  CARD, 
CARDINPUT,  KAL,  KALPAR 
VARI,  KALP 

COMND,  tracking  routine, 
SEARCH 

COMND,  tracking  routine, 

SEARCH 

VARI,  KALP 

VARI,  KALP 

KALPAR,  KAL,  VARI,  KALF, 

CARD,  CARDINPUT 

PAR 3 , KALPAR,  PAR 5 , KAL, 

TIME,  tracking  routines 

COMND,  tracking  routines, 

TIME 

PAR3 , KALPAR,  PAR5 , 

TIME,  tracking  routines 
KAL,  PKQX 

KAL,  PKQX 

VARI,  KALP 

VARI,  CARD,  KALP, 
CARDINPUT,  KAL,  KALPAR 


A 


RMAG 


I 


RMAG 

Distance  from  center  of  earth 
to  the  satellite 

VARI,  KALP,  ORBIT 

RMGDT 

First  derivative  of  RMAG 

VARI,  KALP,  ORBIT 

R0WTK1 

Observer-satellite  vector  in 
topocentric  coordinates 

KALP 

RVDT 

Transverse  component  of  the 
velocity  vector 

VARI,  KALP,  ORBIT 

RX,  RY,  RZ 

ECI  coordinates  of  the 
satellite 

TRANS,  VARI,  KALP,  ORBIT 
PAR3 , PAR5 , KALPAR 

SATNO 

Satellite  number 

CARD,  CARDINPUT,  KAL- 
CARD,  CARDINPUT 

SINI 

Sin  (inclination) 

VARI,  KALP,  ORBIT 

SINLAT 

Sin  (latitude) 

KAL,  KALP,  TRANS,  PAR 3, 
PAR5 , KALPAR 

SINS 

Sin  (sidereal  time) 

TRANS,  ORBIT,  TIME, 
PAR 3 , PAR5 , KALPAR 

SNODE 

Sin  (right  ascension  of  the 
ascending  node) 

VARI,  KALP,  ORBIT 

ST  I ME 

Sidereal  time 

KAL,  ORBIT,  TIME,  PAR3 
PAR5 , KALPAR 

TEMP 

Temporary  storage 

VARI,  KALP 

TEMPL 

Temporary  storage 

VARI,  KALP 

TIME 

Time  into  search 

SEARCH 

TIME2 

Time  of  last  output  to  mount 

MOUNT,  main  tracking 
routines 

TRUEU 

True  anomaly 

VARI,  KALP,  ORBIT 

TSINCE 

Time  since  epoch 

VARI,  TIME,  KALP,  COMND , 
KAL 

VXY 

Constants  for  radius  of 
earth  to  the  equator 

PAR3 , KAL,  PAR5 , KALPAR 

vz 

Constants  for  radius  of 
earth  to  pole 

PAR 5 , KAL,  PAR3 , KALPAR 

XBARK 

Predicted  state  vector  of 
orbital  parameters 

KALP 

XHATK 

Estimated  state  vector  of 
orbital  parameters 

KAL 

YR 

Year  of  card  parameters 

CARD,  COMND,  CARDINPUT, 
KAL,  KALPAR,  KALP 
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SECTION  7 


INTRODUCTION  TO  TV/NOVA  INTERFACE 


The  Kalman  filter  technique  for  updating  orbital  parameters,  described 
in  Section  5 of  this  report,  requires  measurements  of  the  actual  position 
of  the  satellite  in  terms  of  its  azimuth  and  elevation  angles.  A TV  monitor 
on  the  mount  console  displays  the  same  patch  of  sky  at  which  the  telescope 
is  pointed.  The  required  azimuth  and  elevation  measurements  can  be  obtained 
by  means  of  adjustable  cross  hairs  on  the  TV  screen.  This  section  describes 
the  hardware  designed  and  fabricated  to  obtain  the  adjustable  cross  hairs 
on  the  screen  and  to  transfer  the  azimuth  elevation  angular  measurements  to 
the  Active  Tracking  Program  at  operator  command. 

The  camera  interface  positions  two  cross  hairs  on  a video  signal  and 
relays  the  position  of  these  cross  hairs  to  the  CPU  with  an  accuracy  equal 
to  the  resolution  of  the  TV  system. 

To  do  this,  the  interface  must  be  able  to  synchronize  with  the  video 
signal  and  accurately  maintain  its  position  within  the  video  signal. 
Referring  to  Figure  7-1,  one  can  see  the  raster  of  a TV  picture.  We  will 
call  the  solid  raster,  "frame  one"  of  the  picture,  and  dotted  raster, 

"frame  two."  The  two  rasters  are  interlaced  such  that  every  other  line  of 
the  scan  is  done  by  one  frame.  The  slope  of  the  scan  is  such  that  after 
one  complete  horizontal  scan,  the  beam  position  has  moved  down  two  lines 
from  the  top. 

The  interlacing  of  the  frames  is  done  to  provide  as  much  detail  in  the 
picture  as  possible  with  a narrow  bandwidth  signal  and  also  to  avoid  flicker 
on  the  screen. 

'This  interlacing  presents  a problem  to  the  interface.  To  achieve  one 
line  of  resolution,  the  interface  must  draw  half  of  a given  line  on  one 
scan  and  the  rest  of  the  line  on  the  next  scan. 

This  is  pointed  out  in  that,  after  one  half  of  a horizontal  scan,  the 
vertical  position  is  one  line  below.  To  achieve  some  vertical  height  the 
line  must  now  wait  for  one  scan,  start  at  one  line  above  the  end  of  the 
last  line,  and  be  drawn  from  the  halfway  point  to  the  end. 

In  Figure  7-2,  line  #9  is  the  vertical  line  to  be  drawn.  To  do  this, 
the  first  half  of  the  fifth  scan  of  the  frame  #2  raster  is  drawn.  Then  the 
second  half  of  the  fifth  scan  of  the  frame  #2  is  drawn.  In  Figure  7-3,  the 
sixth  line  is  to  be  drawn.  To  do  this,  the  first  half  of  the  fourth  scan 
of  the  second  frame  is  drawn  and  the  second  half  of  the  third  scan  of  the 
first  frame  is  drawn. 


To  accomplish  this  in  hardware,  each  scan  is  divided  into  525  half 
lines.  Each  half  line  is  equal  to  one  vertical  increment.  For  each  horizon- 
tal half  line,  the  vertical  counter  is  incremented  one  count.  On  the 
cross  hairs  can  be  summed  as  either  a positive  (black)  cross  hair  or  a 
negative  (white)  cross  hair. 

The  digital  position  of  the  joystick  is  fed  to  the  computer  as  azimuth 
and  elevation  angles  with  some  scaling  factor.  The  computer  accepts  the 
data  after  the  operator  pushes  the  transmit  button. 

There  is  only  a simplified  explanation  of  the  circuitry.  It  is  only 
intended  to  give  the  reader  a general  knowledge  of  the  function  of  the 
interface. 
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METRIC  SYSTEM 


BASE  UNITS 


Quantity 

Unit 

SI  Symbol 

Formula 

length 

metre 

m 

... 

■Oft 

kilogram 

k« 

tine 

second 

a 

electric  current 

ampere 

A 

thermodynamic  tomporature 

kelvin 

K 

■mount  of  rubeUnce 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

rad 

... 

•olid  angle 

ateradian 

•r 

... 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m It 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration )/s 

angular  acceleration 

radian  per  second  squared 

nd/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

deniity 

kilogram  per  cubic  metre 

kg'm 

electric  capacitance 

farad 

F 

A-s/V 

electrical  conductance 

siemens 

S 

AfV 

electric  field  strength 

volt  per  metre 

... 

V/m 

electric  inductance 

henry 

H 

V-s/A 

electric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

V/A 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

J 

N-m 

entropy 

joule  per  kelvin 

... 

1* 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

(cycle)/* 

illuminance 

lux 

lx 

Im/m 

luminance 

candela  per  square  metre 

... 

cd/m 

luminous  flux 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

Wb 

A/m 

magnetic  flux 

weber 

V-s 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

i/s 

power 

watt 

W 

pressure 

pascal 

Pa 

N'm 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  ateradian 

Wisr 

specific  heat 

joule  per  kilogram-kelvin 

I'kg-K 

stress 

pascal 

Pa 

Nm 

thermal  conductivity 

watt  per  metre-kelvin 

W/m-K 

velocity 

metre  per  second 

m/s 

viscosity,  dynamic 

pascal-second 

Pas 

viscosity,  kinematic 

square  metre  per  second 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(wive)/m 

work 

joule 

1 

N-m 

SI  PREFIXES 


Multiplication  Factors 

Prefix 

SI  Syr 

i ooo  ooo  ooo  ooo  - in11 

ters 

T 

1 000  000  OOC  « 10' 

gigs 

<; 

M 

1 ooo  ooo  *=  io- 

mega 

1 000  - 10' 

kilo 

k 

100  • 101 

hecto* 

h 

10  - 10' 

delu- 

da 

0.1  - 10-’ 

ded* 

d 

001  - 10-> 

centl* 

c 

0.001  « 10- ’ 

milli 

m 

0 000  001  • 10-* 

micro 

M 

0 000  000  001  • 10" 

nano 

n 

0.000  000  000  001  • 10- » 

pico 

rernln 

P 

0.000  000  000  ooo  001  • 10-" 

f 

ooooooooooooooooooi  «=  in-  •• 

atlu 

a 

' To  be  avoidad  whore  possible 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  In  command,  control,  and  communications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications/  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


